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SUMMARY 


An  investigation  was  underwken  to  determine  the  extent  to 
which  presently  available  analytical  methods  can  be  used  to  evaluate 
rotor  blade  air  loads  and  helicopter  rotor  dynamic  response.  The 
Cornell  Aeronautical  Laboratory  (CAL)  variable  inflow  analysis  de¬ 
veloped  under  Contract  DA  44-177-TC-198  was  combined  with  the 
Sikorsky  blade  aeroelastic  analysis.  Appropriate  trim  conditions  were 
developed  and  used  to  apply  the  combing  analysis  to  the  H-34  and 
HU  -  lA  aircraft.  Extensive  test  data  have  been  obtained  for  these  air¬ 
craft.  These  data  were  used  for  a  correlation  study  to  evaluate  the 
accuracy  of  the  analysis.  Flight  conditions  considered  were  steady- 
state  level  flight  at  airspeeds  up  to  113  knots. 

This  correlation  study  showed  that  the  combined  analysis  gave 
significant  improvement  in  accuracy  over  previous  analytical  tech¬ 
niques  using  constant  inflow.  The  definition  of  higher  harmonic  par¬ 
ticipation,  notably  the  blade  root -shear  forces  making  up  the  airframe 
excitation,  was  substantially  better  than  that  obtained  using  constant 
inflow.  In  general,  the  improvement  in  correlation  of  air  loads  and 
moments  was  greater  at  low  airspeeds  than  at  higher  airspeeds  where 
variable  inflow  has  less  effect.  The  variable  inflow  analysis  did  not 
produce  peak-to-peak  stress  values  matching  the  good  correlation  ob¬ 
tained  using  a  constant  inflow  analysis  at  the  high  speed  conditions 
used  for  blade  design.  Although  fairly  good  correlation  with  test 
data  was  obtained,  further  improvement  can  be  made.  The  most  sig¬ 
nificant  area  for  improvement  appears  to  be  in  air-load  definition, 
particularly  in  inflow  definition.  A  further  refined  inflow  analysis  com¬ 
bined  with  the  Sikorsky  aeroelastic  analysis  should  provide  improved 
air  loads  and  moment  values. 

The  analysis  proved  to  be  readily  applicable  to  both  the  HU-IA 
and  H-34  aircraft.  In  addition,  the  analysis  is  finding  application  in 
the  design  of  future  aircraft  now  under  consideration. 
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INTRODUCTION 


There  is  a  real  need  in  the  helicopter  industry  for  an 
aeroelastic  analysis  capable  of  predicting  accurately  the  dynamic 
behavior  of  a  rotor  system.  Such  an  analysis  is  necessary  to 
give  the  industry  a  reliable  method  with  which  to  evaluate  and 
refine  blade  designs.  This  research  program  for  a  ’’Coordin¬ 
ated  Aeroelastic  Research  Program  for  Helicopter  Rotor  Systems” 
was  undertaken  for  the  United  States  Army  Aviation  Materiel 
Laboratories  under  Contract  No.  DA  44-177-AMC-87(T)  to  de¬ 
termine  to  what  extent  the  dynamic  behavior  of  a  rotor  can  be 
predicted  by  application  of  analytical  methods  presently  avail¬ 
able.  A  variable  inflow  analysis,  developed  by  Cornell  Aero¬ 
nautical  Laboratories  (CAL)  to  give  better  definition  of  rotor 
inflow,  was  combined  with  the  Sikorsky  Aircraft  aeroelastic 
analysis,  suitably  modified  to  accept  variable  inflow  and  to  in¬ 
corporate  additional  improvements.  This  analysis  is  a  fully 
coupled  flatwise-edgewise- torsional  analysis  based  on  an  exten¬ 
sion  of  Myklestad’s  method  for  rotating  beams.  The  availabil¬ 
ity  of  both  air -load  and  blade  moment  data  from  flight  tests  of 
both  the  H-34  and  HU-IA  aircraft  (References  I,  6,  9,  and  10) 
provided  a  basis  for  correlation  of  the  analytical  results.  Since 
these  two  aircraft  have  significantly  different  rotor  systems, 
both  aircraft  were  considered  in  order  to  demonstrate  the  abil¬ 
ity  of  the  analysis  to  evaluate  varied  types  of  rotors  and  to  ob¬ 
tain  comparative  data  for  correlation  of  different  rotors. 

This  report  contains  the  results  of  this  program.  The 
validity  of  the  modified  Myklestad  analysis,  which  forms  the 
basis  of  the  Sikorsky  aeroelastic  analysis  is  first  substantiated. 
The  procedures  required  to  carry  out  this  aeroelastic  analysis, 
the  method  of  incorporating  the  variable  inflow  analysis,  and 
the  modifications  made  to  the  analysis  are  then  explained.  Re¬ 
sults  are  given  for  a  study  of  flight  test  data  to  determine  what 
parameters  are  important  in  the  analysis,  and  to  determine  what 
flight  conditions  should  be  used  as  a  basis  for  evaluating  the 
analysis.  Comparisons  between  test  data  and  the  analysis  for 
these  cases  are  presented.  Several  additional  modifications  are 
made  to  the  analysis  in  an  attempt  to  provide  improved  air- load 
definition  through  the  introduction  of  additional  flexibility 
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effects  in  the  calculations.  The  effects  of  these  modifications  are 
presented.  These  results  are  analyzed  to  evaluate  the  ability  of  the 
analytical  procedure  to  predict  dynamic  response  and  to  provide  re¬ 
commendations  for  further  improvements. 
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SUBSTANTIATION  OF  METHODS 


A.  BACKGROUND  OF  ANALYSIS 

Results  of  earlier  correlation  studies  for  the  Sikorsky 
blade  aeroelastic  analysis  have  been  discussed  in  References 
5,  11,  and  12.  Comparisons  of  analysis  (using  constant  inflow) 
with  flight  test  data  indicated  that  calculated  performance  and 
one-half  peak-to-peak  flatwise  blade  stresses  gave  good  agree¬ 
ment  with  flight  measured  values  at  higher  airspeeds.  Satis¬ 
factory  correlation  was  achieved  for  both  articulated  rotor  heli¬ 
copters  (Sikorsky  H-34  and  S-61)  and  rigid-rotor  helicopters 
(Lockheed  CL-475).  However,  it  was  also  found  that  the  assunp- 
tion  of  constant  induced  velocity  imposed  serious  limitations 
on  the  analysis.  Reference  4  concluded  that  better  definition 
of  induced  velocity  was  required  in  order  to:  (1)  predict  higher 
harmonic  air  loads  and  blade  stresses  with  greater  accuracy 
at  all  airspeeds;  (2)  correlate  one-half  peak-to-peak  blade 
stresses  below  100  knots;  and  (3)  determine  analytically  inflight 
fuselage  vibration  levels  at  all  airspeeds. 

Under  this  program,  CAL’s  variable  inflow  analysis 
(Reference  3)  has  been  combined  with  the  Sikorsky  aeroelastic 
analysis  (References  5  and  12)  to  provide  better  definition  of 
induced  velocity  and  associated  higher  harmonic  air  loads.  The 
question  that  remained  for  substantiation  was  whether  the  method 
of  blade  analysis  would  now  adequately  define  higher  harmonic 
blade  bending  moments  and  shears  corresponding  to  these  higher 
harmonic  air  loads.  To  provide  this  substantiation,  Sikorsky 
Aircraft  proposed  to  demonstrate  adequacy  of  the  present  blade 
analysis  method  by  comparing  responses  of  the  ’uncoupled  non- 
rotating  analysis  with  deflections,  slopes,  moments,  and  shears 
as  determined  by  an  independent  closed-form  analytic  solution 
for  sinusoidal  inputs.  Selection  of  a  nonrotating  beam  analysis 
for  substantiation  was  based  upon  the  following;: 

1.  There  was  no  knovvn  closed-form  exact  solution  for 
predicting  the  vibratory  response  of  a  rotating  beam. 

2.  A  nonrotating  beam  comparison  would  place  greater 
demands  on  the  Sikorsky  blade  analysis  method,  since 
higher  harmonic  forced  response  at  a  particular 
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frequency  requires  participation  of  a  greater  number  of 
modes  than  for  the  corresponding  rotating  beam  case. 

The  problems  of  predicting  the  higher  frequency  response 
of  a  1  otating  beam  were  recognized  early  by  Sikorsky  Aircraft  in 
the  development  of  the  Sikorsky  blade  aeroelastic  analysis.  The 
decision  was  made  to  use  a  method  based  upon  direct  inverse  solu¬ 
tion  of  the  blade  dynamic  equations,  rather  than  a  normal  mode  type 
solution.  This  provides  the  advantage  of  full  participation  of  the 
higher  modes(to  as  many  degrees -of -freedom  as  considered),  in 
contrast  to  a  method  which  considers  only  limited  mode  participa¬ 
tion.  The  normal  mode  approach  does  offer  certain  advantages  of 
simplification  of  the  mathematical  formulation  which  can  result  in 
savings  in  computer  time.  Ho»vever,  for  accuracy  the  direct  solu¬ 
tion  approach,  used  by  Sikorsky  Aircraft  in  this  study,  should  equal 
or  surpass  the  normal  mode  approach. 


B.  ANALYTICAL  COMPARISON  OF  METHODS 

A  test  problem  was  set  up  for  determining  the  adequacy  of 
the  present  Sikorsky  blade  analysis  method  for  predicting  higher 
harmonic  blade  response.  A  uniform  pinned-free  beam  with  mass, 
stiffness  and  length  properties  which  approximate  those  of  the  H-34 
main  rotor  blade  was  chosen.  The  beam  is  subjected  to  uniformly 
distributed  oscillatory  loadings  at  frequencies  of  one  per  rev  [203 
cycles  per  minute  (cpm)l  ,  six  per  rev  (1218  cpm),  and  eleven  per 
rev  (2233  cpm).  Objective  of  the  analysis  is  to  compare  responses 
of  the  uncoupled  nonrotating  Sikorsky  blade  analysis  with  corres¬ 
ponding  deflections,  slopes,  moments,  and  shears  as  determined  by 
an  independent  closed-form  analytic  solution. 

The  independent  closed-form  analytic  solution.  Appendix  I, 
was  derived  from  the  classical  Bernoulli -Euler  nonhomogeneous 
differential  equation  for  transverse  vibration  of  a  nonrotating  beam. 
Solution  of  the  equation  gave  resulting  steady- state  harmonic  de¬ 
flections.  slopes,  moments,  and  shears  for  the  uniformly  loaded 
pinned-free  beam. 
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Independent  closed -form  equations  (45).  (46),  (47),  and  (48) 
were  next  programmed  for  the  IBM  7090  computer.  Solutions  were 
obtained  for  p  =  203,  1218,  and  2233  cpm. 

For  comparison,  the  Sikorsky  blade  aeroelastic  analysis  was 
run  in  the  computer  with  identical  input  under  corresponding  nonro¬ 
tating  conditions.  Here, ft  (rotor  angular  velocity)  wes  made  equal 
to  zero,  as  was  blade  twist  in  order  to  decouple  flatwise  and  edgewise 
motions.  In  contrast  to  a  closed-form  solution,  the  blade  for  the 
aeroelastic  analysis  may  not  be  treated  as  a  continuous  structure  but 
must  be  subdivided  into  discrete  segments.  Twenty-four  blade  seg¬ 
ments  were  selected. 

Results  of  the  comparison  are  presented  in  Figures  1  and  2, 
which  give  plots  of  dynamic  response  of  the  blade  as  determined  by 
both  analyses.  Shown  are  resulting  deflections,  moments,  and  shears 
for  sinusoidal  unit- load  excitation  at  frequencies  of  203,  and  2233  cpm  re¬ 
spectively.  The  first  five  modes  for  the  nonrotating  beam  have  na¬ 
tural  frequencies  of  153,  496,  1040,  1765,  and  2700  cpm.  These  are 
helpful  in  the  following  discussion  and  interpretation  of  results. 

Response  of  the  beam  to  1/rev.  excitation  (203  cpm)  is  shown 
in  Figure  1.  This  frequency  is  about  50  cpm  above  the  first  bending 
mode,  about  300  cpm  telow  the  second  mode.  As  expected,  the  de¬ 
flection  curve  is  first  mode  in  character.  Close  agreement  is  noted 
between  analyses  for  resulting  steady-state  harmonic  response  of 
deflections,  moments,  and  shears. 

High  frequency  or  eleven  per  rev  (2233  cpm)  vibratory  response 
is  plotted  in  Figure  2.  In  this  case  frequency  of  response  is  bounded 
by  the  fourth  mode  at  1765  cpm  and  the  fifth  mode  at  2700  cpm.  The 
resulting  deflection  response  curve  exhibits  some  degree  of  fifth- 
modal  response.  For  this  case,  good  agreement  was  also  noted  be¬ 
tween  the  two  analyses,  not  only  for  deflections  but  for  higher  deriva¬ 
tive  bending  moments  and  shears  as  well. 

Similar  correlation  was  found  for  an  excitation  frequency 
of  1218  cpm.  The  results  show  that  the  analysis  is  capable  of 


EXACT  SOLUTION 
SIKORSKY  AIJALYSIS 


100  200  300 
BLADE  RADIUS,  IN. 


100  200  300 
BLADE  RADIUS,  IN. 


100  200  300 
BLADE  RADIUS,  IN. 

FIGURE  1  .  COKPARISON  OF  SIKORSKY  BLADE  ANALYSIS 
WITH  EXACT  SOLUTION  FOR  A  NOinOTATL’IG  EL/vM  • 
CYCLIC  TIP.  MIIIUTL  EXCITATION  F  U  .ILXY  . 


—  EXACT  SOLUnON  I 
— ‘•SIKORSKT  ANALTSI5  “ 


BLADE  RADIUL,  IN. 

FIGURE  2.  CONPARISON  OF  SIKGRSKT  BUDE  ANALYSIS 
VQTH  EXACT  SOLUTION  FOR  A  IJOinCTATHiG  BIAK. 

:  233  CYCLKS  nr.  KiriUTE  dxitaticn  Fn^uzucT . 


7 


giving  an  accurate  definition  of  blade  dynamic  behavior. 


ESTABLISHMENT  OF  a^OTOR  PARAMETERS  AND 
FUSELAGE  TRIM  CONDITIONS 


A  study  was  made  of  the  available  data  to  determine  what 
parameters  affected  the  dynamic  response  most  significantly, 
how  the  program  should  be  modified  to  incorporate  these  para¬ 
meters,  and  exactly  what  flight  conditions  should  be  evaluated 
to  determine  the  validity  of  the  assumptions  made.  Results 
from  extensive  use  of  the  Sikorsky  aeroelastic  analysis  as  well 
as  flight  test  data  were  used  to  determine  what  parameters  and 
modifications  should  be  considered. 

The  flight  test  data  available  for  the  H-34  and  HU-IA 
aircraft  (References  1,  6,  9,  and  10)  were  an  important  source 
of  information  for  the  study.  These  data  served  to  establish  the 
importance  of  a  number  of  parameters,  and  also  as  the  standard 
with  which  analytical  results  could  be  compared.  Perhaps  in 
this  regard,  too  much  may  be  expected  of  the  data.  The  general 
signature  of  the  variation  of  the  air  load,  moment,  and  other 
measurements  with  azimuth  are  probably  given  accurately  by 
the  test  data.  However,  the  test  data  used  were  not  corrected 
for  dynamic  effects.  The  resolution  and  accuracy  of  the  in 
strumentation  and  the  data  recording,  reading,  and  reduction 
methods  used  may  introduce  some  error.  Therefore,  exact 
correlation  with  test  data  is  not  necessarily  the  desired  goal. 
However,  any  substantial  deviation  between  analysis  and  test 
would  indicate  inadequacies  in  the  analytical  procedure. 

There  were  several  areas  in  which  significant  moaifi- 
cations  to  the  Sikorsky  analysis  (References  5  and  12)  were  re¬ 
quired.  These  modifications  include  a  more  realistic  method 
for  obtaining  the  hub  forces  and  momen  .  required  for  trim  and 
more  accurate  damping  coefficients  to  reflect  nonlinear  damping 
rates.  These  are  described  elsewhere  in  this  report. 

The  original  method  used  for  trimming  the  rotor  to  give 
the  proper  hub  moments  and  force  involved  the  use  of  thrust 
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moments  which  are  the  moments  of  the  blade  thrust  about  the  flapping 
axis.  The  steady  and  first  harmonic  components  must  balance  the 
steady  hub  force  and  steady  hub  moment.  This  procedure  locates  the 
center  of  pressure  on  the  rotor  system  correctly,  but  does  not  en¬ 
sure  that  the  values  of  shear  at  the  root  are  those  required  to  balance 
the  hub  forces  and  moments.  In  the  dynamic  system,  it  is  these 
shear  forces  which  must  balance  the  hub  forces.  For  this  reason, 
the  program  was  modified  to  provide  for  trimming  the  rotor  system 
for  the  proper  shear  values. 

Improvement  in  the  method  of  determining  lag-damper  co¬ 
efficients  in  the  analysis  was  also  required.  Examination  of  test 
data  for  lag  motion  at  the  hi'*)  and  the  corresponding  moments  im¬ 
posed  by  the  lag  damper  revealed  that  the  damping  rate  was  a  non¬ 
linear  function  of  rate  of  change  of  lag  angle.  The  exact  manner  in 
which  appropriate  linearized  coefficients  were  introduced  for  each 
response  harmonic  is  described  in  PROVISIONS  FOR  LAG  DAMPER 
on  page  26 

Similarly,  the  lift,  drag,  and  pitching  moment  curves  lead  to 
aerodynamic  damping  te;  ms  which  are  also  nonlinear.  Suitable 
values  for  these  coefficients  were  calculated  as  described  in  AERO¬ 
DYNAMIC  DAMPING  on  page  21  .  Accuracy  in  the  choice  of  aero¬ 
dynamic  damping  coefficients  is  most  important  in  this  analysis. 

TTiey  serve  the  purpose  of  correcting  the  air  loads  on  the  blade  for 
the  effect  of  blade  response.  The  damping  coefficients  are  defined 
by  the  rate  of  change  of  air  load  with  velocity.  These  damping  co¬ 
efficients  must  be  generated  for  use  in  the  forced  response  portion 
of  the  analysis  since  the  aerodynamic  section  does  not  account  for 
blade  flapping  or  flexibility,  but  calculates  air  loads  based  on  a  rigid 
blade  at  the  steady  coning  angle. 

The  data  for  the  H-34  aircraft  showed  that  there  were  large 
edgewise  moments  at  the  outer  end  of  the  blade.  Preliminary  analyses 
gave  no  indication  that  such  moments  existed.  This  difference  was 
reconciled  by  consideration  of  the  effect  of  the  counterweights  in  the 
H-34  blade  spar.  'Phe  counterweights  are  mounted  in  rubber.  Counter¬ 
weight  centrifugal  forces  are  reacted  in  part  by  the  spar  through  the 
rubber,  and  in  part  by  a  moment  reaction  at  the  tip  cap.  In  the  H-34 
test  blade,  most  of  the  force  apparently  was  reacted  by  the  tip  cap, 
for  inclusion  of  this  moment  reaction  would  account  for  the  difference 
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between  measured  and  computed  moments.  Figures  3  and  4  show 
the  effects  of  inclusion  of  corrected  lag-damper  coefficients  and 
counterweight  moment  on  the  edgewise  bending  moments.  Good 
correlation  is  obtained  outboard.  While  inclusion  of  these  effects 
inboard  does  improve  the  correlation,  further  improvement  woukl 
be  desirable. 

The  above  were  the  important  changes  made  to  the  analysis. 
Certain  additional  changes  in  introducing  flexibility  effects  in  the 
air  load  calculations  were  made  for  the  modified  cases.  These  are 
reported  in  MODIFICATIONS  TO  ANALYSIS  on  page  63  .  Examin¬ 
ation  of  the  H-34  test  data  and  other  knowledge  of  the  behavior  of 
the  rotor  system  led  to  the  choice  of  the  following  conditions  for 
evaluation  of  the  analytical  techniques.  For  the  H-34, flight  condi¬ 
tions  consisted  of  trim  level  flight  with  neutral  C.G.  at  41,  70  and 
112  knots.  An  aft  C.G.  condition  was  evaluated  at  73  knots,  while 
the  Bell  HU-IA  rotor  was  evaluated  at  113  knots.  A  speed  of  41 
knots  was  chosen  because  the  bending  stresses  on  the  rotor  blades 
appear  to  reach  a  maximum  point  near  that  speed;  and  this  case 
should  further  reveal  the  characteristics  of  the  analysis  at  low 
speeds.  One  hundred  and  twelve  knots,  on  the  other  hand,  represent 
a  high  speed  case  where  variable  inflow  might  not  be  expected  to 
yield  substantial  improvement  over  a  constant  inflow  approach. 
Seventy  knots  represents  an  intermediate  airspeed  for  which  vari¬ 
able  inflow  is  important.  At  this  speed  the  effects  of  variation  in 
C.G.  should  be  evident,  if  they  are  in  fact  substantial.  The  H-34 
test  data  indicate  that  C.G.  effects  on  air  loads  and  moments  are 
not  large.  As  the  results  will  indicate,  the  analysis  shows  the  same 
lack  of  sensitivity  to  C.  G. 


DEVELOPMENT  OF  AEROELASTIC  MODEL 


The  mathematical  model  was  developed  for  the  study  by  ex¬ 
tending  the  Sikorsky  Aircraft  blade  aeroelastic  program  (References 
5  and  12).  The  basic  analysis  is  fully  described  in  these  references 
and  only  the  significant  extensions  to  the  analysis  will  be  covered  in 
detail  here.  These  extensions  include:  (1)  incorporation  of  the 
variable  induced  velocity  method  of  Reference  3,  (2)  inclusion  of 
rotor  trim  criteria  which  provide  for  hub  pitching  and  rolling  mo¬ 
ments  as  rotor  head  boundary  conditions.  (3)  a  more  accurate  treat¬ 
ment  of  aerodynamic  damping  due  to  blade  flapping,  (4)  provision  for 
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nonlinear  lag-damper  characieristics,  (5)  consideration  of  effects 
of  flexible  blade  motions  on  aerodynamic  load  calculations. 

Before  proceeding  to  describe  the  various  elements  of  the 
analysis  in  detail,  a  general  description  of  the  steps  required  to  ob¬ 
tain  a  solution  v\ill  be  given.  Tlie  procedure  is  outlined  in  Figure 
5.  Block  numbers  used  in  this  section  refer  to  that  figure. 

For  the  helicopter  rotor  to  be  analyzed,  blade  geometry, 
mass,  and  structural  properties  are  required.  In  addition,  fuselage 
drag,  rotor  angular  speed,  aircraft  velocity,  gross  weight,  fuselage 
lift,  and  auxiliary  propulsion  must  be  supplied,  indicated  by  Block 
l.as  well  as  blade  steady-state  two-dimensional  airfoil  data. 

Following  the  solid  lines  in  Figure  5,  initial  rotor  trim  con¬ 
ditions  are  computed,  (Block  2).  Correct  inclination  and  collective 
pitch  of  the  rotor  disk  must  be  determined  to  provide  enough  pro¬ 
pulsive  force  to  overcome  the  net  drag  of  the  aircraft  and  to  support 
the  aircraft.  There  must  also  be  sufficient  cyclic  pitch  to  keep  the 
rotor  in  equilibrium. 

These  initial  trim  conditions,  which  may  be  obtained  by 
applying  a  constant  inflow  condition  to  the  aeroelastic  analysis,  are 
used  in  the  variable  inflow  analysis  (Block  3)  to  determine  the  var¬ 
iable  induced  velocity  distribution  on  the  rotor.  A  discussion  of 
the  variable  inflow  analysis  is  given  in  the  CORNELL  VARIABLE 
INFLOW  PROGRAM  on  page  27  . 

For  determination  of  aerodynamic  loads  (Block  4),  two  di¬ 
mensional  airfoil  data  are  used.  Compressibility  effects  are  taken 
into  account  by  using  Cl  versus  a  ,  Cq  versus  a  .  and  versus 
a  curves  for  angles  of  attack  from  zero  to  thirty  degrees,  to  define 
a  family  of  Cl,  Cq,  and  C|y^  curves  for  Mach  numbers  up  to  .  95.  To 
define  stall  regions  above  angles  of  attack  of  thirty  degrees,  single 
Cl  versus  a  ,  C[^  versus  a,  and  Cm  versus  a  curves  are  taken 
from  Reference  4.  The  blade  to  be  considered  is  subdivided  into 
21  segments.  For  each  of  24  fifteen -degree  azimuth  intervals,  the 
blade  is  analyzed  at  two  blade  pitch  angles.  These  angles  bracket 
the  anticipated  blade  pitch  angles.  Blade  segment  aerodynamic  lifts 
are  then  computed,  from  which  the  moment  of  the  thrust  about  the 
flapping  hinge  is  calculated  as  a  function  of  blade  pitch  angle  and 
azimuth.  The  rotor  blade  pitch  necessary  to  maintain  the  rotor 
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system  in  equilibrium  is  calculated  by  a  thrust  moment  iteration 
technique,  described  in  PROVISION  FOR  ROTOR  TRIM  on  page  17, 
which  forces  the  first  harmonic  thrust  moment  about  the  flapping 
hinge  to  converge  to  the  required  value  (Block  5).  These  values  for 
blade  pitch  angles  are  obtained  by  a  linear  interpolation  procedure. 
Blade  segment  aerodynamic  lift,  drag,  and  pitching  moments  are 
recalculated,  using  the  previously  trimmed  values  for  pitch,  at  each 
of  the  21  radial  blade  locations  as  harmonic  functions  of  azimuth 
angle  (Block  6).  The  magnitude  of  calculated  helicopter  gross 
weight  is  compared  with  the  required  gross  weight,  where  calculated 
helicopter  gross  weight  is  based  upon  the  mean  value  of  thrust  mo¬ 
ment.  rotor  lift,  rotor  drag,  and  the  blade  radial  location  of  the 
resultant  thrust  vector.  If  the  two  differ  by  more  than  the  gross 
weight  tolerance  (normally  1%),  refinements  to  rotor  lift  and  rotor 
drag  and  other  pertinent  rotor  attitude  parameters  are  made  (Block 
7).  Previously  determined  blade  pitch  angles  are  incremented  and 
resubstituted  to  initiate  new  rotor  equilibrium  calculations.  TTie 
iteration  procedure  is  continued  until  the  specified  tolerance  on 
gross  weight  is  satisfied.  Flatwise  aerodynamic  damping  coeffi¬ 
cients  are  then  calculated  (Block  8). 

An  additional  refinement  to  rotor  trim  has  been  incorporated 
which  employs  steady  hub  rolling  and  pitching  moments  as  boundary 
conditions.  An  analytical  discussion  of  this  trim  procedure  is 
given  in  PROVISION  FOR  ROTOR  TRIM  on  page  17  .  The  dashed 
lines  in  Figure  5  trace  the  procedure.  Expressions  for  the  hub 
rolling  and  pitching  moments  can  be  written  in  terms  of  blade  root 
shears  and  moments;  thus  flexible  dynamics  as  well  as  aerodynamic 
effects  are  accounted  for  in  the  trim  procedure.  For  an  articulated 
rotor,  steady  and  first  harmonic  flatwise  root  shear  and  steady 
edgewise  root  shear  arc  employed. 

The  steady  edgewise  shear  can  be  determined  from  rotor 
torque  requirements.  To  determine  the  flatwise  root  shears  an  un¬ 
coupled  flatwise  Myklestad  analysis  is  used  (Block  9).  This  gives 
the  approximate  shear  values  for  the  given  aerodynamic  loading. 

This  analysis  is  used  instead  of  a  fully  coupled  analysis  to  decrease 
the  running  time  of  the  program  for  the  shear  trim.  A  3  by  3  matrix 
is  formulated  which  relates  increments  of  collective  and  cyclic 
blade  pitch  to  increments  of  steady  and  first  harmonic  blade  root 
shear.  The  increments  of  collective  and  cyclic  pitch  obtained  from 
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the  solution  of  the  matrix  equation  are  the  values  needed  to  satisfy 
the  specified  trim  conditions  (Block  10). 

Calculation  of  the  aerodynamic  damping  coefficients  used  in 
the  Myklestad  equations  to  account  for  blade  flapping  is  discussed 
in  full  in  AERODYNAMIC  DAMPING  on  page  21  .  Sine  and  cosine  N 
per  revolution  flatwise  perturbation  velocities  are  superimposed 
upon  the  rotor  equilibrium  condition  determined  by  the  thrust  mo¬ 
ment  iteration  scheme.  Changes  in  the  harmonics  of  radial  distri¬ 
bution  of  blade  segment  lift,  drag,  and  pitching  moment  form  three 
N  by  N  complex  aerodynamic  damping  matrices,  where  N  indicates 
the  number  of  harmonics  of  rotor  speed  considered.  At  present, 
the  main  diagonal  elements  of  the  three  N  by  N  matrices  are  intro¬ 
duced  into  the  fully  coupled  Myklestad  aeroelastic  blade  program. 
Interharmonic  coupled  damping  coefficients  are  not  considered. 
Complex  aerodynamic  damping  coefficients  are  then  multiplied  by 
the  appropriate  blade  segment  velocities  within  the  Myklestad  analy¬ 
sis.  This  procedure  thus  defines  the  harmonics  of  complex  aerody¬ 
namic  damping  force.  Provisions  are  made  in  the  Myklestad  analy¬ 
sis  to  include  aerodynamic  damping  forces  because  blade  flapping 
velocities  are  not  included  in  defining  an  inflow  angle  for  air- load 
calculations.  The  pitch  angle  refinement  calculation  from  the  flat¬ 
wise  root -shear  matrix  method  does  include  the  effect  of  the  damping 
forces  and  thus  ensures  that  further  air- load  calculations  do  include 
the  influence  of  blade  segment  flapping  velocities. 

Convergence  of  the  shear  iteration  technique  is  obtained  when 
a  prescribed  root-shear  tolerance  condition  is  met.  After  calcula¬ 
tion  of  the  edgewise  and  torsional  damping  coefficients  (Block  11), 
the  harmonics  of  air  loads  associated  with  the  trim  condition  are 
applied  to  the  fully  coupled  Myklestad  relations.  This  analysis  is 
described  in  detail  in  Reference  12.  The  lag-damper  coefficients 
used  in  this  analysis  are  supplied  as  input  data  and  are  obtained  by 
the  method  described  in  PROVISIONS  FOR  LAG  DAMPER  on  page  26. 

In  the  present  study,  steady  boundary  conditions  were  used 
in  the  Myklestad  analysis.  However,  it  is  possible  to  incorporate 
the  airframe  dynamic  response  in  the  analysis  and  use  matched 
rotor  and  airframe  impedance  at  the  hub  as  boundary  relations.  The 
dynamic  interaction  of  rotor  system  and  airframe  then  would  be  con¬ 
sidered  directly  in  the  analysis.  The  basis  for  this  approach  has 
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been  developed  by  Sikorsky  Aircraft  (Reference  5). 

The  procedure  to  this  point  is  that  used  for  the  basic  calcu- 
laiion  for  the  five  flight  conditions  evaluated.  Modifications  were 
made  to  this  procedure  in  an  attempt  to  improve  air -load  definition. 
The  changes  involve  a  repetition  of  the  entire  procedure,  as  indi¬ 
cated  by  the  circled  solid  lines,  recalculating  the  air  loads  to 
account  for  the  effects  of  blade  bending  from  the  previous  solution. 
Blade  segment  torsional  rotation  and  flatwise  slopes  are  introduced 
into  the  air- load  calculations. 

The  variable  induced  velocity  distribution  used  in  the  aero¬ 
dynamic  analysis  can  be  refined  by  introducing  flexible  motions 
from  the  fully  coupled  Myklestad  analysis  into  the  Cornell  analysis. 

These  modifications  are  described  in  MODIFICATIONS  TO  ANALYSIS 
on  page  63 . 

The  procedure  followed  after  computation  of  the  new  air 
loads  is  exactly  the  same.  The  rotor  system  is  retrimmed  using 
the  new  air- load  distribution.  The  fully  coupled  Myklestad  analysis 
is  then  used  to  determine  new  values  for  blade  response.  Additional 
iterations  may  be  made  if  desired  for  further  refinement. 


A.  PROVISION  FOR  ROTOR  TRIM 

A  thrust  moment  iteration  scheme  has  been  used  in  the  Si¬ 
korsky  aeroelastic  analysis  for  trimming  a  rotor  system  to  obtain 
collective  and  cyclic  blade  pitch  for  a  given  steady- state  level  flight 
condition.  Hiis  approach  equates  the  steady  thrust  and  moments 
produced  by  the  air- load  distribution  over  the  rotor  with  the  re¬ 
quired  thrust  and  moment  on  the  rotor  hub.  Since  the  forces  and 
moments  transmitted  to  the  rotor  hub  depend  directly  upon  the  shear 
and  moments  at  the  inboard  end  of  the  rotor  blades,  a  new  shear 
trim  iteration  scheme  was  developed  to  refine  the  pitch  values  given 
by  the  thrust  moment  iteration  scheme.  The  shear  iteration  incre¬ 
ments  blade  pitch  to  produce  the  required  shear  at  the  inboard  end 
of  the  blade. 

Initial  rotor  trim  conditions  are  established  by  nneans  of  the 
thrust  moment  iteration  procedure,  in  which  rotor  equilibrium  is 
deter.mined  by  enforcing  sufficient  blade  pitch  to  give  required  first 
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harmonic  thrust  moment  about  the  flapping  hinge,  and  setting  the 
steady  thrust  moment  equal  to  the  following  value. 

TMq  =  x(gw  -  Lfus  *  Lrotor  )  (Steady  thrust  moment) 

TcoinaTccw  TZ  (1) 

where 

X  is  the  radial  location  of  thrust  vector 

gw  is  the  gross  weight  of  ship 

Lfyg  is  the  lift  of  fuselage 

Lrotor  is  vertical  component  of  rotor  drag 

b  is  the  number  of  blades 

a  (  is  the  angle  of  tip  path  plane 

Po  is  the  blade  rigid  body  coning  angle 

As  a  result  a  3  by  3  symmetric  thrust  moment* blade  pitch 
matrix  is  established  which  allows  direct  solution  of  the  blade  pitch 
angle.  Development  of  the  matrix  is  based  upon  defining  the  Fourier 
coefficients  of  thrust  moments  as  functions  of  linearized  relations 
between  thrust  moment  and  blade  pitch.  Final  converged  values  of 
blade  pitch  are  obtained  by  specifying  a  tolerance  on  the  calculated 
aircraft  gross  weight.  Azimuthal  thrust  moment  calculations  using 
blade  segment  theory  are  made  for  an  assumed  upper  and  lower 
blade  pitch  setting.  A  linear  variation  of  thrust  moment  and  blade 
pitch  for  various  azimuth  positions  is  assumed. 

TM(«^)  =  m(«^)  [Aq+Bi  sin  ^  +  A  j  cos  s'  ]  '♦'  b(  i//  )  (2) 

where 

TMliJ;)  is  the  azimuthal  thrust  moment 

m(^)  is  the  slope  of  azimuth  thrust  moment  blade  pitch 
lines 
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sin  ^  +A1CO6  ^  is  the  blade  pitch  angle 

b(^  )  is  the  intercept  of  the  thrust  moment  blade  pitch 

zero  line 


The  azimuthal  distribution  of  thrust  moment  may  be  written 
in  the  following  harmonic  form. 

TM(^)  =  TMq  +  TMi  cosf  sin^  +  . . .  (3) 

where 


TMq  is  the  steady  thrust  moment 

TMj  is  the  first  harmonic  cosine  component  of  thiust 
moment 

TWi  is  the  first  harmonic  sine  component  of  thrust 
moment. 


Upon  comparing  the  coefficients  of  a  Fourier  series  expres¬ 
sion  of  Equation  (2)  with  Equation  (3)  the  following  3  by  3  matrix 
results. 


2TMo 


1  df 


.a. 


®11  1  ®12  I  “13 

“21 J  *22  j  “23 


Bl 


(4) 


TMj 

[™1  |_a3i  j  «32  j  “33j 

where  the  matrix  coefficients  are  defined  as; 

aji  812  *  “21  df 

“22  >  cos^f  df  aj3  -  831  « I  Aifi(f  )8inf  df 

“33  sin^fdf  833  *  832  <^K(f  )8inf  coef  df 


Inversion  of  the  matrix  will  yield  the  collective  (Aq)  and 
cyclic  pitch  (A^.  B^)  for  rotor  trim.  Tie  to  the  nonlinear  charac¬ 
teristics  of  thrust  moment  variation  with  blade  pitch,  several  iter¬ 
ations  on  blade  pitch  must  be  made  to  meet  gross  weight  and  flapping 


4 

1 

f 
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requirements. 


Once  convergence  is  obtained,  the  shear  trim  criteria  are 
applied  to  obtain  refined  values  for  blade  pitch.  A  set  of  boundary 
relations  associating  fixed  hub  rolling  and  pitching  moments  and 
vertical  force  as  functions  of  blade  root  steady  and  first  harmonic 
flatwise  and  edgewise  shears,  moments,  and  torques  is  used  to 
determine  the  shear  required  at  the  blade  root.  Boundary  relations 
may  be  chosen  to  allow  the  use  of  any  type  of  blade  root  construction, 
whether  fully  articulated,  semi -articulated,  or  rigid,  and  also  allow 
a  rotor  system  to  be  trimmed  for  off  neutral  center  of  gravity  con¬ 
dition.  For  a  fully  articulated  rotor,  the  following  trim  identities 
are  used  to  define  the  shears  required  at  the  blade  root. 

F,  =  Sn^b  cos  Or  (steady  vertical 

^  hub  force)  (6) 

Mxx  =  sin  QiSq  b  +  e  cos  Qj'S’i  b  (steady  hub  rolling 

— 5 —  moment)  (7) 


(steady  hub  pitching 
moment)  (8) 


Vertical  hub  force  is  given  by  the  rotor  lift  required.  Steady 
hub  rolling  and  pitching  moments  may  be  obtained  from  an  aero¬ 
dynamic  equilibrium  analysis.  With  an  estimated  value  for  shaft 
angle,  at.  the  required  steady  edgewise  shear,  So^.  and  steady 
and  first  harmonic  flatwise  shear,  So^,  Sl^,  and^l^,  are  deter¬ 
mined.  Actual  Sq^  on  the  blade  is  determined  from  rotor  torque 
requirements  and  actual  flatwise  shears  Sq.  Sj,  and  Si,  are  found 
from  applications  of  given  air  loads  to  an  uncoupled  flatwise  Mykle- 
stad  blade  analysis.  Changes  in  blade  pitch  required  to  satisfy  the 
shear  trim  conditions  can  then  be  determined  with  an  iterative  pro¬ 
cedure  using  the  following  equation. 


Fs  '"-s] 
^0  ^0 

^Sq  I  I 

^0 

aAq  I 

— 

Sl^-Si 

1  1 

Al 

5^-5, 

AAq  I  ^A 1  1  AB  1 

aSj  1  j 

Bl 

i 

_AAo  ’  AAl  !  ABi  _ 
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The  difference  elements  of  the  matrix  are  obtained  by  sep¬ 
arately  incrementing  collective  pitch  and  sine  and  cosine  values  of 
cyclic  pitch,  and  determining  the  corresponding  flatwise  root  shears 
from  the  flatwise  blade  analysis. 

The  incremental  values  of  collective  and  cosine  and  sine 
components  of  blade  pitch  are  aAqi  AAi,  and  aBi* 

Upon  inverting  the  matrix  new  values  of  blade  pitch  are  ob¬ 
tained  from 

0  )  =  (Ao AAo)  +  (Aj^  AAi)  cos^  +  (Bi  +  aBi)  sin^ 

(10) 

where  Aq,  Ai,  Bi,  are  blade  pitch  angles  established  by  the  thrust 
moment  iteration  scheme  or  from  earlier  passes  of  the  shear  iter¬ 
ation. 


b.  aerodynamic  damping 

A  method  of  analysis  was  developed  for  determining  aero¬ 
dynamic  damping  coefficients  associated  with  rotor  blade  flapf^ng 
motions.  The  method  makes  use  of  a  perturbation  technique  where¬ 
in  increments  are  made  to  the  blade  segment  flapping  velocities 
associated  with  a  particular  rotor  trim  condition.  From  this  pro-  , 
cedure,  corresponding  changes  in  air  loads  are  obtained.  The 
aerodynamic  damping  coefficients  acting  on  an  element  of  the  blade 
can  then  be  defined  by  terms  of  the  form 


AL  ^ 

’  l-^\  • 

/  apm 

p  1 

‘  l~®p  /  ‘ 

w 

These  coefficients  represent  the  respective  change  in  blade  lift, 
drag,  and  pitching  moment  for  unit  change  in  blade  flapping  velocity. 

The  total  perturbation  velocity  is  introduced  in  the  aero- 
elastic  model  by  incrementing  the  normal  component  of  the  total 
blade  segment  velocity  by  the  following  n/rev  harmonic  series: 
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(11) 


(  AUp)i  =  A(r  4  )i  = 

where 

(  AUp)i  =  A(r  4  )i 


!(»■  4  )ni  +  (r  4  )ni 

is  the  increment  of  normal  compo¬ 
nent  of  total  blade  segment  velocity 


i  is  the  blade  segment  location 


4 


is  the  flapping  angular  velocity 
is  the  blade  segment  radial  location 


n  is  the  number  of  harmonics 

s 

(r  /3  )jjj  is  the  harmonic  of  sine  perturbation 
velocity 

.  Q 

(r  0  is  the  harmonic  of  cosine  perturba¬ 

tion  velocity 


Damping  coefficients  for  each  harmonic  of  the  fully  coupled 
Myklestad  blade  response  analysis  must  be  obtained.  Incrementing 
(Up)j  by  each  harmonic  of  flapping  velocity  as  expressed  in  Equation 
(if)  yields  the  following  complex  matrix  of  aerodynamic  damping 
coefficients: 


where  A  is  a  general  notation  indicating  air  load  (lift,  drag,  or 
pitching  moment) 


is  the  complex  aerodynamic  damping  coefficient 
matrix 

is  the  real  part  of  complex  aerodynamic  damping 
coefficient  matrix 


AA  ~]1  is  the  imaginary  part  of  complex  aerodynamic 
Ar  4.  Ji  damping  coefficient  matrix 
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Real  and  imaginary  parts  of  the 


ping  coefficient  matrix  are  expressed  as: 

R 

[  AA  1  .  1 
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complex  aerodynamic  dam- 
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(13) 


where 
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(14) 

is  the  incremental  change  in  the  nth/rev  cosine 
component  of  total  air  load  due  to  a  nih/rev  cosine 
component  of  flapping  velocity 

is  the  incremental  change  in  the  nth/rev  sine  com¬ 
ponent  of  total  blade  segment  air  load  due  to  the 
nth /rev  cosine  component  of  segment  flapping  vel¬ 
ocity. 


^(r  0)li 


AA 


s 

ni 


A(r  0  ) 


ni 


is  the  incremental  change  in  the  nth/rev  cosine  com¬ 
ponent  of  total  blade  segment  air  load  due  to  nih/rev 
sine  component  of  segment  flapping  velocity 

is  the  incremental  change  in  the  nih/rev  sine  com¬ 
ponent  of  total  blade  segment  air  load  due  to  nth /rev 
sine  component  of  segment  flapping  velocity 


In  the  fully  coupled  Myklestad  analysis,  only  the  diagonal 
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elements  in  the  coefficient  matrix  have  been  included.  Interharmonic 
damping  coefficients,  which  are  defined  by  the  off  diagonal  elements 
have  been  neglected.  These  coefficients  were  neglected  on  the  basis 
that  their  effects  would  be  small.  However,  examination  of  later  re¬ 
sults  indicated  that  the  effect  of  interharmonic  damping  becomes  im¬ 
portant  when  large  harmonic  blade  flapping  occurs.  A  further  dis¬ 
cussion  is  given  in  CORRELATION  OF  PREDICTED  AND  FLIGHT 
TEST  LOADS  on  page  33. 

It  was  noted  that  the  coefficients  on  the  main  diagonal  were 
almost  identical  in  magnitude.  Thus  the  first  harmonic  damping 
coefficient  was  substituted  for  the  remaining  main  diagonal  coef¬ 
ficients.  Similar  identities  were  also  noted  for  the  off  diagonal  co¬ 
efficients. 


The  elements  in  the  complex  aerodynamic  damping  matrix 
obey  the  law'  of  complex  differentiation  of  an  analytic  function.  If  it 
is  assumed  that  the  Cauchy -Riemann  relationships  apply,  the  following 
must  hold: 
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Calculation  of  each  element  did  reveal  that: 


(15) 


aa 


ni 


AA® 

III 


A(r  0  A(r  0  )nj 

Because  the  terms  were  not  exactly  equal,  an  average  of  the 
two  values  was  used.  For  simplification,  define  the  elements  of  the 
complex  aerodynamic  coefficient  matrix  by: 
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are  indices  indicating  change  in  nth  harmonic  of  airload 
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due  to  qth  harmonic  of  flapping  velocity, 
is  the  blade  radial  station. 


Likewise  each  element  of  the  imaginary  part  of  the  coefficient 
matrix  will  be  designated  by: 
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Equation  (15)  will  then  be  expressed  as: 

(Cnq)i  ■  (Cnq)|^  "*■  j  (Cnq)j^  (18) 

The  complex  damping  force  Fj,  will  include  the  influence  of 
all  interharmonic  coefficients.  Its  magnitude  is  obtained  by  multi¬ 
plying  Equation  (18)  by  the  appropriate  harmonic  of  blade  segment  vel¬ 
ocity.  To  define  the  blade  velocity,  suppose  the  displacement  of  the 
blade  is  expressed  by  the  following  harmonic  series: 

<yq)i  =  <yq)r  +J(yq)i'  (19) 

The  velocity  is  obtained  by  taking  the  time  derivative  of 
Equation  (19). 


(Vq)i  =  (Vq)f  +  J(Vq){ 


(20) 


The  complex  aerodynamic  damping  force  is  thus: 

„  <  Fn)i  =  t  (Cnq)i  (Vqlj 

(Fn)i=  I  f  Cnq)i(Vq)r  -  (Cnq),'  (Vq){  1  +J  f(C„q)J(v  )f+ 

0=1  I*-  J  L 


(C„q). 
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force  is;  r  i  1  r  d  o  i  i 

(Fi)f=  [(Cii)f  (vj)F-  (Cii)|  (v,){J+  •  |(Cin),  (Vn)r  -  (Ci„)i(v„)}] 


As  an  example  the  real  part  of  the  first  harmonic  damping(22) 
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where  n  is  the  number  of  harmonics  used. 

When  inierharmonic  damping  coefficients  are  neglected  the 
equation  for  reduces  to 

(Pn)i  "  (^nn^i  ^'^nh 

This  simplification  was  used  in  the  program.  It  was  believed  that 
inclusion  oi  interharmonic  terms  would  significantly  increase  the 
running  time  of  the  program  without  providing  much  of  a  change  in 
the  damping  force.  However,  as  was  stated  earlier,  results  given 
in  CORRELATION  OF  PREDICTED  AND  FLIGHT  TEST  LOADS  on 
page  33,  indicate  that  the  interharmonic  terms  may  be  important 
for  some  cases. 

C.  PROVISIONS  FOR  LAG  DAMPER 


To  incorporate  the  nonlinear  characteristics  of  a  lag-hinge 
damper  with  relief  valve  provision  in  the  fully  coupled  analysis,  the 
follow'ing  procedure  was  developed.  For  a  representative  speed  dur¬ 
ing  trimmed  level  flight,  measured  edgewise  bending  moments  at  the 
blade  root  and  lead- lag  angles  are  used  to  define  the  magnitude  of 
the  damping  coefficients.  Moment  and  lag  angle  data  are  analyzed  to 
obtain  harmonics  in  terms  of  azimuth  angle.  The  harmonics  of  angu¬ 
lar  blade  lagging  velocity  are  then  calculated  from  the  lag  angle  har¬ 
monics.  These  may  be  described  by: 


M^(  «(/  ) 

C  O') 
t  ( i') 


cos  n  ^  +  Mp  sin  n  ^ 
cos  ny>  +  sin  ni|> 
X'oftEnSinn  ^  n  cos  n  ^ 


measured  blade  root 
edgewise  moment  (24) 
measured  blade  lagging 
angle  (25) 

computed  lagging  angu¬ 
lar  velocity  (26) 


From  these  relations,  the  lag  damper  coefficient  is  defined 
in  general  as: 


C( 


(27) 


In  view  of  the  harmonic  definition  of  the  blade  lagging  angular 
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velocity  and  its  root  edgewise  moment,  a  damping  coefficient  is  con¬ 
structed  for  each  harmonic  of  rotor  speed. 


This  coefficient  is  defined  by  the  sine  component  of  root  edge¬ 
wise  moment  divided  by  the  cosine  component  of  lagging  blade  velo¬ 
city  and  the  cosine  component  of  root  edgewise  moment  divided  by 
the  sine  component  of  lagging  blade  velocity. 
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(28) 


The  damping  coefficient  used  to  describe  the  nonlinear  char¬ 
acteristic  of  a  lag-hinge  damper  is  expressed  for  each  harmonic  of 
rotor  speed  where  n  indicates  the  order  of  the  harmonic. 


The  damping  coefficient  for  each  harmonic  of  rotor  speed  is 
then  introduced  in  the  fully  coupled  Sikorsky  blade  response  program 
where  it  is  coupled  to  the  appropriate  harmonic  of  lag  motion. 

n.  THE  VARIABLE  INFLOW  PR(X3RAM 


Description  of  Analysis 


The  variable  inflow  distributions  required  for  the  Sikorsky 
Aeroelastic  Program  were  calculated  using  the  Cornell  Aeronautical 
Laboratory  Variable  Inflow  Computer  Program  developed  under  TRECOM 
(now  USAAVLABS)  Contract  No.  DA  44-177-TC-698  (References 
3,  7.  and  8).  The  mathematical  model  consists  of  the  representatiun 
of  each  blade  by  a  segmented  lifting  line,  and  the  helical  wake  of  the 
rotor  by  a  mesh  of  discrete  segmented  vortex  filaments  consisting  of 
both  trailing  and  shed  vorticity.  The  shed  vortex  elements  result 
from  the  azimuthwise  variation  of  bound  circulation  and  the  trailing 
vortex  elements  result  from  the  spanwise  variation  of  bound  circula¬ 
tion.  The  circulation  of  the  wake  for  each  blade  changes  with  azi¬ 
muth  position  and  is  periodic  for  each  rotor  revolution.  The  blades 
are  divided  into  a  finite  number  of  radial  segments,  and  the  induced 
velocity  at  the  center  of  each  selected  blade  segment  is  computed  by 
summing  the  contributions  of  each  bound  trailing  wake,  and  shed  wake 
vortex  filament.  The  contribution  of  each  vortex  filament  is  obtained 
through  use  of  the  Biot-Savart  equation  which  expresses  the  induced 
velocity  in  terms  of  the  circulation  strength  of  the  vortex  filament 
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and  its  geometrical  position  relative  to  the  blade  segment  at  which 
the  induced  velocity  is  desired. 

The  program  can  be  considered  as  having  two  basic  pans: 

(1)  the  determination  of  the  blade -wake  geometry  and  (2)  the  deter¬ 
mination  of  the  circulation  distribution  and  its  associated  induced 
velocity  distribution.  The  wake  geometry  is  determined  by  assuming 
that  each  element  of  wake  vorticity  travels  with  the  speed  and  in  the 
direction  of  the  instantaneous  velocity  imparted  to  it  at  the  time  it 
was  generated  at  the  blade.  The  wake  geometry  is  thus  fixed  by  the 
choice  of  the  input  wake  velocity  distribution,  rotor  dimensions,  and 
flight  condition.  In  reality,  the  wake  transport  velocities  and  resul¬ 
tant  wake  geometry  are  greatly  influenced  by  interaction  effects  be¬ 
tween  the  various  elements  of  vorticity  in  the  wake.  However,  the 
inclusion  of  such  wake  vortex-vortex  interaction  effects  was  beyond 
the  scope  of  the  variable  inflow  program  used  for  the  results  repor¬ 
ted  herein.  The  hound  circulation  distribution  is  determined  by  re¬ 
lating  the  wake  circulations  to  the  bound  circ  lations,  expressing 
the  wake  induced  velocities  in  terms  of  the  unknown  bound  vortex 
strengths  by  means  of  the  Biot-Savari  equation,  and  developing  a 
set  of  simultaneous  equations  relating  the  bound  circulation  and  local 
blade  angle  of  attack  at  each  blade  segment.  These  equations  thus 
involve  the  knossn  flight  condition,  wake  geometry,  lift-curve  slope, 
and  blade  motion  and  control  parameters  and  the  unknown  bound 
circulation  values.  Solution  of  these  equations  yields  the  desired 
bound  circulation  values  which,  when  combined  with  the  appropriate 
geometrical  relations  in  the  Biot-Savart  law,  produce  the  required 
induced  velocity  distribution. 


Assumptions  of  the  Analysis 

The  state  of  the  art  concerning  variable  inflow  theory  is  con¬ 
tinuously  advancing  and  many  of  the  assumptions  listed  below  are 
currently  under  investigation.  In  this  research  program,  the  Cor¬ 
nell  variable  inflow  analysis  of  Reference  3  was  used.  There  is  a 
newer  version  of  the  Cornell  program  in  which  there  are  provisions 
for  rolling  up  the  trailing  wake  filaments  and  the  elimination  of  the 
shed  wake  filaments  in  the  far  wake  region  (generally  assumed  to 
lie  a  quarter  of  a  revolution  downstream  of  the  blade).  This  treat¬ 
ment  of  the  vortex  interaction  and  dissipation  phenomenon  was  com- 
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pleted  by  the  Cornell  personnel  during  the  period  of  this  3ntract 
but  was  not  included  in  the  analysis  reported  herein.  Since  further 
study  of  the  rotor  wake  geometry  problem  is  currently  in  progress 
at  Cornell  Aeronautical  Laboratory,  and  the  inclusion  of  the  above- 
mentioned  initial  wake  modifications  served  principally  to  reduce 
computer  time  rather  than  to  effect  any  major  changes  in  the  pre¬ 
dicted  induced  velocity  and  air  load  distribution,  the  original  var¬ 
iable  inflow  program  described  in  Reference  3  was  used  in  obtaining 
the  results  presented  herein.  This  method  incorporates  a  wake 
model  for  which  the  near  wake  geometrical  distribution  of  trailing 
and  shed  wake  elements  is  retained  in  the  far  wake.  This  original 
program  had  been  previously  evaluated  by  United  Aircraft  Corpora¬ 
tion  Research  Laboratory  and  Sikorsky  Aircraft,  and  much  of  the 
discussion  and  conclusions  concerning  the  program  are  based  on 
that  evaluation. 

The  following  is  a  list  of  the  principal  assumptions  of  the 
Cornell  variable  inflow  program  used  to  compute  the  induced  velo¬ 
city  distributions  employed  in  this  report. 

1.  The  flight  condition  is  steady. 

2.  The  blades  and  rotor  wake  are  represented  by  vortex 
elements  based  on  lifting  line  theory  rather  than  finite 
vortex  core  and  lifting  su:*face  theory.  Thus,  aerody¬ 
namic  boundary  conditions  at  the  blade  are  satisfied  at 
only  one  point  on  the  chord. 

3.  Viscous  dissipation  effects  on  the  wake  circulation 
strengths  are  neglected  inasmuch  as  the  circulation 
strength  of  a  given  wake  element  is  constant  with  time. 
However,  the  number  of  wake  revolutions  retained  in  the 
analysis  can  be  limited. 

4.  Unsteady  aerodynamic  effects  resulting  from  the  varia¬ 
tion  of  bound  circulation  strength  with  blade  azimuth  po¬ 
sition  are  represented  by  the  inclusion  of  shed  wake  ele¬ 
ments  and  the  provision  for  varying  circulation  strength 
between  successive  trailing  wake  elements.  However, 
due  to  the  rather  coarse  spacing  of  the  shed  vorticity  and 
the  representation  of  the  blade  by  a  lifting  line,  the  ac- 
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curacy  of  the  unsteady  aerodynamic  model  is  questionable. 

5.  Stall  is  included  by  limiting  the  maximum  lift  coefficient 
value  and  thus  the  maximum  circulation.  It  is  assumed 
that  the  stalled  lift  is  circulatory. 

6.  Reversed  flow  effects  are  included  by  the  correct  direc¬ 
tions  of  the  velocity  components.  The  location  of  the 
bound  vortex  does  not  change  from  the  1/4-  to  the  3/4- 
chord  position. 

7.  Small  angle  assumptions  are  included  in  the  analysis. 

8.  Mach  number  and  Reynolds  number  effects  are  included  by 
using  as  input  to  the  program  lift-curve  slope  values 
'Ahich  may  vary  over  the  rotor  disc  but  not  with  angle  of 

niiack. 

9.  In -plane  components  of  induced  velocity  are  neglected. 

10.  Interference  effects  of  the  rotor  hub»  fuselage,  etc.  are 
neglected. 

11.  I’hc  wake  geometry  is  specified  by  wake  transport  veloci¬ 
ties  which  are  assumed  to  be  equal  to  the  velocities  im¬ 
parted  to  the  various  wake  elements  at  the  time  they  were 
generated  at  the  blade.  Wake  vortex  interaction  and  vis¬ 
cous  dissipation  effects  are  thus  neglected. 

12.  Blade  flapping  harmonics  above  the  first,  lag,  and  bending 
effects,  with  the  exception  of  steady  bending,  are  not  in¬ 
cluded  in  the  establishment  of  the  wake  geometry,  and 
therefore  in  the  determination  of  the  blade  element  to 
wake  element  distances. 


E.  COUPLING  OF  THE  VARIABLE  INFLOW  PROGRAM 
VVITH  T(^E  SIKORSKY  AEROELASTlC  ANALVS(S 


The  variable  inflow  program,  in  its  present  form  .separates 
the  air- load  problem  from  the  blade  response  problem.  The  circula- 
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tion  and  induced  velocity  distributions  are  determined  from  known 
blade  motions  and  control  parameters.  If  the  circulation  distribu¬ 
tion  thus  computed  is  used  directly  to  obtain  air  loads,  it  is  general¬ 
ly  found  that  the  aircraft  trim  conditions  (lift,  drag.  etc. )  and  the 
blade  boundary  conditions  are  not  satisfied.  Thus,  in  coupling  the 
variable  inflow  program  to  the  Sikorsky  aeroelastic  analysis,  it  was 
necessary  to  establish  an  approach  which  would  satisfy  these  con¬ 
ditions. 


The  original  approach  considered  involved  an  iteration  pro¬ 
cedure  which  used  only  that  portion  of  the  variable  inflow  program 
defining  tl  3  wake  geometric  coefficients  (  o-’s)  by  the  Biot-Savart 
equation.  An  initial  estimate  of  blade  motions,  control  parameters, 
and  lift  distribution  was  obtained  from  the  blade  dynamic  response 
program  using  constant  momentum  inflow,  but  s'itisfying  the  blade 
boundary  conditions,  and  aircraft  trim  conditions.  The  circulation 
values  were  then  obtained  from  the  lift  distribution  and  combined 
with  the  wake  geometric  coefficients  determined  by  the  variable  in¬ 
flow  program  to  produce  an  initial  estimate  of  the  induced  velocity 
distribution.  This  variable  induced  velocity  distribution  was  then 
used  as  input  to  the  blade  response  program  and  a  new  circulation 
distribution  was  obt  ained.  This  procedure  was  continued  until  the 
circulation  values  converged.  Due  to  the  change  in  inflow  distribu¬ 
tion.  it  was  then  necessary  to  include  an  iteration  to  adjust  the  col¬ 
lective  pitch,  etc.  to  satisfy  aircraft  trim  requirements. 

This  approach  was  discountinued  due  to  convergence  problems 
associated  with  the  circuit  ion  iteration  method  which  were  encoun¬ 
tered  whenever  wake  vortex  elements  were  in  close  proximity  to 
the  blade  stations  at  which  the  air  loads  were  being  evaluated.  The 
nature  of  the  difficulties  can  be  illustrated  by  expressing  the  circula¬ 
tion  iteration  procedure  in  terms  of  a  representative  set  of  simultan¬ 
eous  equations  which  have  the  form  indicated  in  Equation  (6)  of 
Reference  3. 

Tji  =  i  +  ^ki  *  (29) 

This  relation  expresses  the  bexmd*  circulation  values  of  the  jth  itera¬ 
tion.  r  |[.  in  terms  of  the  local  geometric  angle -of- attack  distribution 
and  blade  plunging  motions.  I|c.  the  wake  geometric  coefficients.  ^|ci, 
the  wake  circulation  values  of  the  j-1  iteration,  t  j'l.  and  the  pro- 
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ducts  of  the  lift-cuive  slope  and  blade  semichord,  a^bj^.  Nondimen- 
sionalizing  the  above  equation  by  dividing  by  ,  relating  the 

wake  circulation  values  to  the  bound  circulation  values,  and  writing 
it  in  matrix  form; 


{rj'  =  [i]  + 

Bound  Blade 

Circulation  Dynamics 


V  {r]’’'  (30) 

Wake  Wake  Circulation 

Geometry  in  Terms  of 

Bound  Circulation 


Induced  Velocity 


The  convergence  of  the  iteration  method  represented  by  the 
above  set  of  simultaneous  algebraic  equations  was  investigated  by 
examining  the  effect  of  the  magnitude  of  the  y  matrix  coefficients 
on  convergence  characteristics.  It  was  determined  that,  for  a 
specific  value  of  I,  convergence  difficulties  should  be  expected 
whenever  the  magnitude  of  any  of  the  matrix  coefficients  approaches 
approximately  .5  to  1.  These  values  correspond  to  critical  9 
coefficients  of  .5/a|^bj^  to  1/aj^bj^.  Examination  oi  the  elements 
of  the  variable  inhow  program  9  matrices  for  several  flight  con¬ 
ditions  of  interest  indicated  potential  convergence  difficulties  under 
the  following  circumstances; 


1.  Whenever  blade  segment  lengths  were  less  than  approx¬ 
imately  5  per  cent  of  the  rotor  radius. 

2.  At  low  speeds  where  blades  pass  in  close  proximity  to 
vortex  elements  generated  by  previous  blades. 

3.  At  high  speeds  where  reverse  flow  effects  cause  shed 
vortexes  to  approach  and  pass  the  blade  generating  the 
vortexes. 


In  view  of  these  anticipated  difficulties  this  iteration  approach  was 
abandoned. 

The  revised  iteration  method  involves  more  extensive  use  of 
the  existing  variable  inflow  program  in  that  the  procedure  given  by 
the  program  for  calculating  the  rotor  circulation  distribution  is  em- 
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ployed  in  addition  to  the  procedure  for  calculation  of  the  wake  geomet¬ 
ric  coefficients.  In  effect,  this  replaces  the  portion  of  the  original 
circulation  iteration  method  by  a  closed-form  solution,  thereby  eli¬ 
minating  the  Iteration  difficulties  previously  mentioned.  It  is  noted 
that  the  method  for  computing  circulation  includes  the  assumption  of 
constant  lift-curve  slope  values  for  each  Mach  number.  However, 
as  mentioned  previously,  the  Cornell  variable  inflow  solution  does 
not  satisfy  blade  boundary  conditions,  and  hence  an  iteration  between 
the  variable  inflow  program  and  the  blade  analysis  is  necessary.  It 
should  also  be  noted  that  aircraft  trim  conditions  (lift,  drag,  etc. )  are 
satisfied  by  the  blade  analysis  through  adjustments  in  collective  pitch, 
cyclic  pitch,  and  rotor  angle  of  attack. 

Due  to  the  dependence  of  the  variable  induced  velocity  distri¬ 
bution  obtained  from  the  variable  inflow  program  on  the  blade  motion 
and  control  parameter  input,  it  was  necessary  to  include  a  test  of  the 
convergence  of  this  input  after  each  pass  through  the  blade  aeroelastic 
analysis.  However,  it  was  generally  found  that  inflow  convergence 
was  rapid,  and  cne  iteration  pass  was  sufficient. 


CORRELATION  OF  PREDICTED  AND  FLIGHT  TEST  LOADS 

The  analysis  was  applied  to  the  mai:  rotors  of  the  H-34  and 
HU-IA  helicopters.  Four  flight  conditions  were  evaluated  for  the 
H-34,  while  one  flight  condition  was  considered  for  the  HU-IA. 
These  cases  correspond  to  steady  state  level  flight  conditions  for 
which  flight  test  data  are  available.  The  cases  considered  are: 

Case  1.  H-34;  steady-state  level  flight  at  41  knots  with 
neutral  C.  G. 

Case  2.  H-34;  steady-state  level  flight  at  73  knots  with  aft 
C.G. 

Case  3.  H-34;  steady- state  level  flight  at  70  knots  with 
neutral  C.  G. 

Case  4.  H-34;  steady-state  level  flight  at  112  knots  with 
neutral  C.  G. 

Case  5.  HU-IA;  steady-state  level  flight  at  113  knots  with 
neutral  C.  G. 

Test  data  for  the  cases  were  (Stained  from  the  following 
sources;  Case  1  from  Table  V  in  Reference  9;  Case  2  from  Table 


33 


104  in  Reference  10;  Case  3  from  Run  2  in  Reference  6;  Case  4 
from  Table  21  in  Reference  10;  and  Case  5  from  condition  31  in 
Reference  1. 

The  following  sections  present  the  input  data  used  in  the 
analysis,  the  significant  results  obtained,  and  the  correlation  with 
test  data. 

A.  BASIC  BLADE  DATA 


The  H-34  and  HU-IA  have  significantly  different  rotor  sys¬ 
tems.  TTie  H-34  has  a  four-bladed  fully  articulated  rotor  system 
with  lag  dampers,  llns  helicopter  was  evaluated  at  a  gross  weight 
of  11,  800  pounds.  The  HU-IA  has  a  two-bladed  teetering  rotor 
system.  Since  the  blade  is  not  free  to  hunt,  lag  dampers  are  not 
required.  This  helicopter  was  evaluated  at  a  gross  weight  of  6027 
pounds. 


The  basic  blade  data  for  the  two  aircraft  are  presented  in 
Tables  1  and  2  (page  35  and  36  ).  These  tables  give  the  radial 
station  breakdown  used  in  the  aeroelastic  analysis,  and  the  blade 
properties  at  each  station.  Ttie  calculated  natural  frequencies  of 
these  blades  are  shown  in  Table  3  (page  37  ). 

Presented  in  Table  4  (page  38  )  are  the  trim  conditions  used 
for  each  of  the  five  cases. 

Trim  values  of  first  harmonic  flatwise  root  shear  for  the 
H-34  cases  were  obtained  from  calculations  of  required  steady  hub 
rolling  and  pitching  moments.  A  preliminary  aerodynamic  analysis 
based  on  constant  inflow  was  used  to  obtain  the  required  rolling 
and  pitching  moments.  Zero  first  harmonic  shear  was  used  for 
the  HU-IA. 

Analysis  of  the  teetering  rotor  of  the  Bell  HU-lA  requires 
the  choice  of  appropriate  boundary  conditions  for  each  harmonic 
considered.  Tlie  boundary  conditions  chosen  correspond  to  those 
given  in  Reference  2.  For  response  to  steady  loads,  the  rotor  was 
considered  cantilevered  both  flatwise  and  edgewise.  For  the  1, 

3,  5,  and  7  per  rev.  harmonics,  the  blade  was  pinned  flatwise  and 
cantilevered  edgewise,  ror  the  2,  4,  and  6  per  rev.  harmonics 
the  blade  was  cantilevered  flatwise  and  pinned  edgewise.  A  canti- 
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TABLE  1 

BASIC  BLADE  DATA  H-34  HELICOPTER 
GROSS  WEIGHT  11,  800  POUNDS.  TWIST  -8  DEGREES 
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TABLE  2 

BASIC  BLADE  DATA  HU-IA  HELICOPTER 
GROSS  WEIGHT  602"  POUNDS.  TWIST  -12  DEGREES 
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10. 175  581 


TABLE  3 


CALCULATED  BLADE  NATURAL  FREQUENCIES 


Natural  Freq^uencies,  CPM 
flatwise  Edgewise  Flatwise  Edgewise 
Mcxie  Pinned  Pinned  Cantilevered  Cantilevered 


HU-IA^ 

1 

574.2 

715.3 

375.5 

530.5 

2 

1063. 1 

1886.6 

1126.8 

3198.  8 

3 

1718.5 

3679. 6 

2237.  9 

H-34** 

1 

872.3 

1857.6 

2 

1607.1 

5616.0 

3 

2819.7 

•  D.  =  311.5  RPM 
*•0.  =  216  RPM 
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BASIC  FLIGHT  PARAMETERS 
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levered  condition  was  assumed  for  torsion. 

The  H-34  was  considered  to  be  pinned  both  flatwise  and  edge¬ 
wise  and  cantilevered  in  torsion. 

Data  for  the  0012  airfoil  of  the  H-34  were  available  at  Sikor¬ 
sky  Aircraft.  Cl,  Cq,  and  Cjyi  versus  ct  curves  were  used  for  the 
range  of  Mach  numbers  up  to  .  95.  Data  tor  the  HU-IA  0015  airfoil 
were  supplied  by  Bell  Helicopter  Company.  Only  Cl  and  Cq  daw 
were  available,  so  0012  values  were  used  for  the  Cw.  The  aifference 
in  Cjyj  values  between  the  0012  and  0015  airfoil  was^lieved  to  be 
relatively  insignificant  in  this  analysis. 


B.  VA.’^IAPLE  INFUDW  ANALYSIS  INPUT 


The  input  parameters  required  for  the  variable  inflow  analysis 
consist  of  rotor  geometry,  flight  conditions,  choice  of  blade  segments, 
wake  geometry,  airfoil  data,  and  blade  motion  and  control  parameters. 

The  rotor  dimensions  and  flight  condition  parameters,  with 
the  exception  of  the  rotor  angle  of  attack,  were  obtained  from  Refer¬ 
ences  1,  6,  9,  and  10. 

The  lift-curve  slope  values  were  obtained  by  averaging  the 
slope  values  in  the  unstalled  linear  region  of  the  Sikorsky  0012  air¬ 
foil  data  and  the  0015  data  supplied  by  the  Bell  Helicopter  Company 
for  the  appropriate  Mach  numters  at  each  point  in  the  rotor  disc.  An 
angle -of- attack  value  of  11.  5  degrees  was  used  to  limit  the  wing  lift 
coefficient  in  the  stalled  region  for  both  configurations,  the  H-34 
and  the  HU-IA.  A  check  of  the  Mach  number  angle -of -a tuck  combi¬ 
nations  at  each  point  in  the  rotor  disc  revealed  that  these  assumptions 
had  little  effect  on  the  results  because  the  flight  conditions  chosen  for 
this  study  were  essentially  unstalled;  that  is,  all  blade  settions  for 
the  HLI-IA,  with  the  exception  of  those  in  the  proximity  ci  the  rever¬ 
sed  flow  region  boundary  at  the  higher  advance  ratios,  operated  in  ‘ 
the  linear  lift-curve  slope  region  and  below  the  assumed  11.5  degree 
stall  angle,  whereas  the  H-34  operated  slightly  into  the  stalled  region 
above  the  11.5  degree  stall  angle. 

In  Reference  7  the  Cornell  variable  inflow  analysis  was  used 
to  calculate  air  loads  for  the  HU-IA  at  110  knots.  Flight  test  data 
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for  Condition  67  of  Reference  1  were  used  for  correlation  (and  not 
Condition  31  which  was  used  in  the  present  study. )  The  results  in¬ 
dicated  a  stalled  region  for  this  condition.  The  absence  of  stall  in 
the  present  study  is  due  to  the  use  of  aircraft  gross  weight  for  rotor 
lift.  Actual  lift  was  somewhat  higher  as  can  be  seen  from  Figure  37 
which  shows  the  steady  air -load  component  for  the  HUIA  at  113 
knots.  Since  rotor  lift  was  underestimated,  low  values  for  pitch 
angle  were  calculated  (see  Table  5  on  page  46  )  which  indicated  that 
the  rotor  did  not  stall.  In  the  work  of  Reference  7,  the  problem  of 
determining  actual  lift  was  avoided  by  using  pitch  angle  and  not  lift 
as  the  input  parameter.  The  effects  of  downwash  on  the  fuselage  and 
negative  fuselage  lift  due  to  forward  till  of  the  aircraft  are  most 
likely  the  reason  for  rotor  lift  exceeding  aircraft  gross  weight. 

Nine  blade  radial  segments  were  used  to  represent  both  the 
H-34  and  the  HU-IA  blades  in  the  variable  inflow  program.  Tlie 
segment  midpoints  were  located  at  r  /R  =  .  23, .  375,  .  525,  .  65,  .  75, 

.  85,  .  925,  .  965,  and  .  99.  Tlie  number  and  p>ositioning  of  the  seg¬ 
ments  determine  not  only  the  blade  segment  distribution  but  also  the 
number  of  wake  segments,  their  location  and  spacing.  Results  com¬ 
paring  the  harmonics  of  blade  loading  for  a  five  and  nine  segment  re¬ 
presentation  of  a  model  blade  at  an  advance  ration  of  .  15  in  Refer¬ 
ence  3  indicate  that  the  choice  of  the  number  of  segments  can  affect 
the  air  loads.  To  evaluate  this  effect,  the  use  of  seven  and  nine  seg¬ 
ments  was  compared  for  the  H-34  rotor  at  the  41-knot  flight  condi¬ 
tion  (Figure  6).  A  comparison  of  the  harmonics  of  induced  velocity 
and  air -load  distribution  for  the  two  segment  distributions  produced 
relatively  small  differences  except  at  the  tip  segment,  and  all  dif¬ 
ferences  were  less  than  the  difference  between  the  analytical  results 
and  the  model  data.  However,  nine  segments  were  used  to  give  a 
reasonable  up  segment  distribution  and  to  allow  for  possible  increased 
sensitivity  at  higher  speeds. 

An  azimuth  increment  of  15  degrees  was  employed  in  defining 
the  rotor  and  wake  model  in  the  variable  inflow  program  and  was  the 
minimum  value  possible  considering  program  limits.  A  comparison 
of  the  harmonics  of  inflow  obtained  from  the  variable  inflow  program 
using  15-degree,  22.  5-degree,  and  30-degree  azimuth  increments 
for  a  39-knot  flight  condition  indicated  that  the  results  for  the  22.5- 
degree  and  15-degree  increments  are  essentially  the  same.  It  was 
not  possible  to  investigate  the  effect  of  small  increments.  However, 
it  is  possible  that  the  treatment  of  the  shed  wakes  in  the  immediate 


40 


7 

9  Si^xa-kJirs 


r  'R  -  .985 


90  180  270 
AZIMJ7H  KMle.,  Di:XRil£S 
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CGRNm.  VAKIaBLS  lifLOW  PROGRAM  ON  BLADS 


vicinity  of  the  blade  could  produce  difficulties  if  small  azimuth  incre¬ 
ments  are  used. 

The  number  of  wake  revolutions  retained  for  the  wake  model 
for  the  H-34  helicopter  41 -knot,  70-knot,  and  112-knot  conditions 
were  three,  two  and  one,  respectively.  Also,  one  wake  revolution 
was  retained  for  the  HU-IA  helicopter  113- knot  condition.  The  choice 
of  number  of  wake  revolutions  was  based  on  the  determination  of  the 
sensitivity  of  the  resultant  air  loads  to  this  parameter  as  indicated  by 
the  results  of  Reference  1.  It  is  mentioned  in  Reference  3  that  trun¬ 
cating  the  wake  at  three  revolutions  provided  sufficient  accuracy  for 
model  rotor  at  an  advance  ratio  of  .  15.  The  use  of  three  revolutions 
was  also  found  to  be  adequate  for  the  H-34  rotor  at  an  advance  ratio 
of  .  11.  The  insensitivity  of  the  air  loads  to  the  effects  of  additional 
wake  revolutions  for  these  flight  conditions  is  due  to  the  increased 
average  distance  between  this  portion  of  the  wake  and  the  rotor.  In¬ 
asmuch  as  the  rotor  is  moving  away  from  its  wake  even  faster  at 
higher  speeds,  it  is  to  be  expected  that  at  the  higher  flight  speeds  the 
number  of  wake  revolutions  necessary  to  maintain  the  same  accuracy 
will  be  decreased.  This  reasoning  was  used  to  determine  the  smaller 
number  of  wakes  used  for  the  H-34  at  70- knot  and  112-knot  conditions, 
and  the  HU-lA  at  113-knot  condition.  The  distance  between  corres¬ 
ponding  points  on  the  rotor  and  the  farthest  wake  revolution  for  these 
conditions  was  maintained  at  approximately  the  same  value  as  that  of 
the  41 -knot,  three  revolutions  case.  The  H-34  at  112  knot  condition 
was  investigated  for  both  one  and  two  wake  revolutions.  Differences 
in  calculated  inflow  were  very  small. 

The  wake  transport  velocity  input  to  the  variable  inflow  pro¬ 
gram  was  generally  assumed  to  be  the  constant  induced  velocity  deter¬ 
mined  from  momentum  considerations.  It  is  realized  that  this  assump¬ 
tion  is  not  completely  accurate  inasmuch  as  significant  distortions  of 
the  helical  wake  structure  can  be  produced  by  wake  vortex  interaction 
effects.  The  use  of  this  model  is  justified  for  the  following  reasons: 

(1)  No  provision  is  made  in  the  existing  variable  inflow  program  for 
putting  in  realistic  time  varying  wake  transport  velocities,  (2)  Ques¬ 
tions  still  exist  as  to  what  displacement  time  history  should  be  used 
for  a  particular  rotor  and  operating  condition  (Reference  3),  and  (3) 
Initial  results  presented  in  Reference  3  also  indicate  that  the  use  of  a 
uniform  inflow  distribution,  determined  from  momentum  considera¬ 
tions,  resulted  in  computed  air  loads  which  correlate  reasonably  well 
with  measured  air  loads  for  a  two-bladed  rotor. 
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Calculated  variable  inflow  distributions  for  the  H-34  at  41  and 
il2  knots  are  shown  in  Figures  7  and  8.  These  plots  show  the  inflow 
distribution  over  the  rotor  disc.  Tliere  is  a  wide  variation  in  inflow 
over  the  disc  with  steep  velocity  gradients  near  the  outer  edge. 


C.  ROTOR  TRIM  ANGLES 


For  the  calculation  of  the  results  presented  here,  it  was  ne¬ 
cessary  to  choose  estimated  values  for  blade  pitch  and  shaft  inclina¬ 
tion  for  the  determination  of  rotor  inflow  in  the  variable  inflow  analy¬ 
sis.  Application  of  this  inflow  to  the  Sikorsky  aeroelastic  analysis 
results  leads  to  a  redefinition  of  these  parameters  for  proper  rotor 
trim.  Close  correlation  between  the  values  used  to  determine  inflow 
and  the  resulting  values  given  by  the  aeroelastic  analysis  was  obtained. 

Listed  in  Table  5  (page  46)  is  a  comparison  between  the  analy¬ 
tical  and  flight  test  values  for  collective  and  cyclic  pitch,  rotor  shaft 
angle,  and  steady  coning.  The  pitch  values  are  given  with  reference 
to  the  root  for  the  H-34  and  to  station  23  for  the  HU-IA.  Positive 
signs  indicate  pitch  leading  edge  up,  shaft  angle  rearward,  and  coning 
up. 


D.  PRESENTATION  OF  RESULTS 


Results  of  the  analysis  have  been  plotted  with  comparable 
flight  test  data.  Plots  representing  the  azimuthal  variation  of  airloads, 
flatwise  bending  moments,  edgewise  bending  moments,  torsional  mo¬ 
ments,  and  control  loads  have  been  made  and  are  available,  but  are 
not  presented  in  this  report.  Plots  are  presented  for  the  radial  varia¬ 
tion  of  air  loads,  flatw'iae  bending  moments,  and  edgewise  bending  mo¬ 
ments.  These  plots  show  the  harmonic  content  of  the  loads.  These 
were  obtained  from  a  24-point  harmonic  analysis  of  the  data.  The  har¬ 
monic  coefficients  are  dehned  by  the  following  relations. 


where 


Y(^)=  1 

Acn  ■  I 
12 


7 

Aco  J] 

n=l 

i  y(*k) 

k=0 


cos  n  A 


(31) 
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COMPARISON  OF  ROTOR  TRIM  ANGLES  FROM  ANALYSIS  AND  TEST 
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sn  = 


T2  2]  YC'I'k)  sin 
K=0 


The  value  of  the  plotted  variable  at  azimuth  angle  is  Y(  ), 
and  the  angle  ^  k  defined  by  ^^=(1510®.  To  aid  in  comparing 
the  harmonic  content  of  measured  and  analytical  results,  the  harmon¬ 
ics  are  tabulated  in  terms  of  amplitude  and  phase  angle  8n.  de¬ 
fined  by 


or 


An  cos  (n  -0n)  *  ^cn 


^sn 

_Asn^' 

^cn 


Agn  sin  ni// 


(32) 


The  following  sign  convention  was  used  in  the  presentation  of 
results.  A  positive  value  of  section  aerodynamic  loading  indicates 
a  force  in  the  lift  direction.  Positive  signs  indicate  compression  in 
the  upper  surface  of  the  blade  for  flatwise  bending  moment,  compres¬ 
sion  in  the  leading  edge  of  the  blade  for  edgewise  bending  moment,  a 
couple  tending  to  rotate  blade  trailing  edge  upward  for  torsional  mo¬ 
ment,  and  an  upload  on  the  end  of  the  pitch  horn  for  pitch  horn  load. 


E.  DATA  CORRELATION 


Examination  of  the  results  shown  in  Figures  19  through  45 
reveals  that  fairly  good  correlation  between  test  and  analysis  was  ob¬ 
tained.  However,  significant  differences  do  exist,  which  indicate 
that  further  work  is  required  to  produce  an  analysis  capable  of  re¬ 
liably  predicting  blade  behavior.  In  general,  there  appear  to  be  no 
significant  differences  between  the  correlation  found  for  air  loads  and 
moments.  Computed  moment  distributions  correspond  to  those  which 
would  be  expected  from  the  applied  air- load  distributions,  which  shows 
the  dynamic  system  is  responding  properly.  The  following  evaluation 
summarizes  the  significant  points  concerning  data  correlation. 


Case  1.  H-34  at  41  Knots 

This  is  a  low  speed  case  where  substantial  discrepancies  had 
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been  found  between  analytical  results  based  on  constant  inflow  and 
test  data.  Use  of  variable  inflow  does  improve  correlation.  Exam¬ 
ination  of  the  azimuthal  variation  in  air  loads  shows  that  the  analytical 
results  do  follow  the  general  signature  of  the  test  data.  However, 
particularly  i^ear  the  blade  tip.  the  rapid  changes  in  air  loads  are  not 
matched  by  analysis.  This  indicates  that  the  higher  harmonic  content 
of  the  analytical  air  loads  is  not  sufficiently  defined.  Examination  of 
the  radial  dis  ribution  of  the  harmonics  (Figures  19  through  21)  re¬ 
veals  that  the  steady  and  first  harmonic  amplitude  and  phase  match 
H-34  test  data  quite  well.  Somewhat  larger  differences  occur  in  the 
second  harmonic  amplitude  outboard  of  r/R  =  .  7.  The  third,  fourth, 
and  fifth  harmonics  show  differences  in  amplitude  and  in  phase  near 
the  tip.  The  present  analysis  does  not  produce  an  adequate  definition 
of  the  air  loads  near  the  tip. 

T^e  radial  flatwise  moment  plots  of  Figures  22  through 

24  bear  this  out.  Low  values  of  a*r  loads  on  the  outboard  half  of 

the  blade  were  calculated  for  the  second,  third,  and  fourth  harmonics. 
It  is  these  harmonics  of  flatwise  moment  which  were  found  to  be  low. 
Belter  correlation  was  found  for  the  first  and  third  harmonics  where 
air  loads  correlate  weil  also.  The  reason  for  the  difference  between 
steady  moments  cannot  be  determined.  Reference  10  indicates  that 
the  steady  flatwise  moments  which  are  tabulated  are  not  reliable, 
since  consistent  values  of  steady  moment  could  not  be  obtained  for 
the  same  flight  condition  on  different  days.  This  was  due  to  difficul¬ 
ties  in  obtaining  a  true  zero  reference.  The  azimuthal  plots  of  flat¬ 
wise  moments  confirm  the  above. 

Chordwise  moment  radial  distribution,  shown  in  Figures 

25  through  27,  follow  the  same  pattern.  Torsional  moments  and  con¬ 
trol  loads  calculations  yielded  higher  values  for  the  lower  harmonics 
than  was  found  in  test. 


Cases  2  and  3.  H-34  at  73  Knots  with  Aft  C.G.  and  at 
70  Knots  with  Neutral  C.G. 

These  two  cases  should  be  considered  together.  They 
differ  mainly  in  the  C.G.  location  which  produces  a  difference  in 
steady  moment  acting  on  the  rotor  head,  and  consequently  different 
amounts  of  first  harmonic  flapping.  Flapping  should  produce  only 
small  changes  in  airloads,  as  shown  by  the  test  data.  However,  sig- 
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nificant  differences  in  air  loads  were  found  in  the  analytical  results. 
Investigation  revealed  that  these  differences  were  largely  due  to  the 
absence  of  interharmonic  damping  in  the  analysis.  This  investiga¬ 
tion  is  discussed  more  fully  later  in  this  section.  The  main  point  is 
that  the  presence  of  first  harmonic  flapping  produces  a  substantial 
difference  in  the  second  harmonic  content  of  the  air  loads.  For  the 
aft  C.  G.  case  the  magnitude  of  the  second  harmonic  is  approximately 
half  the  magnitude  for  the  neutral  C.  G.  case.  The  difference  is  al¬ 
most  equal  to  the  correction  which  inclusion  of  interharmonic  damping 
would  make. 

Examination  of  the  results  shows  that,  in  general,  correlation 
was  better  for  the  70- knot  neutral  C.  G.  case  than  for  the  aft  C.  G. 
case.  There  is  a  trend  throughout  the  results  of  poor  correlation 
near  an  azimuth  angle  of  0  degrees.  This  difference  will,  of  course, 
affect  all  harmonics.  With  the  exception  of  the  second  harmonic 
difference,  air  loads  for  the  two  cases  correlated  with  test  equally 
well.  Flatwise  moments  showed  somewhat  better  correlation  for 
the  neutral  C.  G.  case  than  for  the  aft  C.  G.  case.  There  is  a  sub¬ 
stantial  difference  in  chordwise  moments.  For  the  neutral  C.  G.  case, 
chordwise  moments  correlate  well  outboard  and  differ  mainly  by  a 
shift  in  the  steady  moment  at  inboard  stations.  The  aft  C.G.  case 
shews  a  large  shift  in  first  harmonic  phase  as  well  as  a  difference 
in  steady  moments  inboard,  with  comparable  correlation  outboard. 
Calculated  torsional  moments  and  control  loads  correlated  fairly 
well  with  test  data  for  the  advancing  blade,  but  not  very  closely  for 
the  retreating  blade. 

■ 

Case  4.  H-34  at  112  Knots 

The  results  for  this  case  are  presented  as  radial  distributions 
of  harmonics  as  shown  in  Figures  28  through  36.  Air  loads  for  this 
case  show  fairly  good  correlation  with  test  data  in  terms  of  am¬ 
plitudes.  However,  there  is  a  significant  shift  in  phase  between  test 
and  analytical  results.  This  shift  in  phase  is  also  evident  in  the  flat¬ 
wise  bending  moment  curves.  These  moments  otherwise  show  good 
correlation  except  at  azimuth  angles  in  the  neighborhood  of  zero  de¬ 
grees.  Comparison  of  the  harmonics  of  air- load  and  flatwise  bending 
moment,  both  in  magnitude  and  phase  angle,  shows  that  the  correla¬ 
tion  of  flatwise  moments  follows  ouite  consistently  the  correlation 
of  air  loads.  The  harmonics  for  which  air  loads  correlate  closely 
with  H-34  data  also  show  close  correlation  of  flatwise  moments.  If 
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calculated  air  loads  are  small,  moments  are  small  also.  This  in¬ 
dicates  that  lack  of  correlation  is  largely  due  to  inaccurate  definition 
of  air  loads.  This  t(^ic  is  discussed  in  EVALUATION  OF  ANALYSIS 
on  page  68.  Chordwise  moments  show  a  sigpiificant  difference  in 
steady  moment  inboard.  Calculated  higher  liar  monies,  especially 
the  third  and  fourth,  are  much  lower  than  H-34  data.  The  general 
signature  of  the  chordwise  moments,  however  is  quite  good.  Tor¬ 
sional  moments  and  control  load  correlation  showed  large  differences 
for  the  retreating  blade. 

Case  5.  HU-IA  at  113  Knots 


Results  for  this  case  are  presented  in  Figures  37  through  45 
which  show  radial  distribution  of  the  air -load  and  bending  moment 
harmonics.  Air- load  correlation  is  fairly  good.  There  is  a  phase 
shift  between  test  and  calculated  data.  This  is  most  evident  near 
zero  azimuth  and,  as  in  other  cases,  leads  to  poor  correlation  of 
air  loads  in  that  region.  Unfortunately,  rotor  thrust  was  underesti¬ 
mated  by  about  15%  (see  page  39)  which  yielded  low  calculated  steady 
air  loads,  and  low  values  for  the  secc^  and  third  harmonics.  The 
radial  plots  show  good  correlation  for  the  first  and  fourth  harmonics. 
Flatwise  moment  calculated  values  were  high  for  the  first  and  fourth 
harmonics  but  other  harmonics  correlated  very  well.  Chordwise  mo¬ 
ment  correlation  was  more  erratic.  In  general,  calculated  moment 
values  were  less  than  test  values.  This  is  most  evident  in  the  steady, 
second,  and  fourth  harmonics.  The  differences  in  the  second  and 
fourth  harmonics  are  most  likely  due  to  the  boundary  conditions  used 
for  these  harmonics.  The  rotor  was  considered  to  be  pinned  for  these 
modes  to  reflect  the  flexibility  of  the  rotor  support  and  the  antisym¬ 
metric  character  of  the  response.  The  results  show  that  higher  mo¬ 
ments  should  exist  on  the  inboard  section  of  the  blade.  This  indicates 
that  a  cantilevered  boundary  condition  would  have  been  more  realistic. 
Torsional  moments  and  control  load  have  approximately  the  correct 
magnitude  but  poor  phase  correlatioii. 

F.  COMPARISON  WITH  CONSTANT  INFLOW 

The  use  of  the  variable  inflow  analysis  has  produced  signifi¬ 
cant  improvements  in  correlation  with  test  data.  Typical  results  for 
the  H-34  at  41  knots  and  at  112  knots  are  shown  in  Figures  9 
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through  12.  At  41 -knots,  constant  inflow  (without  a  tip  loss  factor) 
correlates  poorly  with  H-34  data,  while  the  variable  inflow  results 
show  good  correlation.  On  the  other  hand,  at  112  knots,  there  is 
much  less  difference  beirween  correlation  obtained  with  variable  in¬ 
flow  and  with  constant  inflow.  The  reason  for  this  is  that  at  41  knots 
the  induced  velocities  vary  widely  across  the  rotor  disc,  and  constant 
inflow  is  a  poor  approximation.  However,  for  speeds  in  the  range  of 
112  knots,  the  variation  of  inflow  have  much  less  effect  on  air  loads 
and  the  constant  inflow  approxima«^ion  is  more  valid.  Therefore, 
better  correlation  would  be  expected.  There  is  a  difference  in  phase 
between  calculated  air  loads  and  flight  test  data  as  shown  in  Figure 
11.  Similar  phase  differences  were  found  in  the  work  of  Reference 
7.  It  appears,  therefore,  that  the  phase  difference  results  from  the 
inflow  distribution  used  in  the  analysis. 


C.  ROOT'SHEAR  FORCES 


The  use  of  variable  inflow  has  produced  more  realistic  values 
for  root-shear  forces.  These  are  the  forces  which  are  transmitted 
to  the  rotor  hub  to  produce  airframe  excitation.  Earlier  studies 
(Reference  5)  showed  that  root  shear  forces  predicted  by  constant 
inflow  were  less  than  those  required  for  airframe  excitation.  Calcu¬ 
lated  airframe  response  using  these  shear  forces  was  much  below 
measured  response. 

Plots  of  the  harmonics  of  root- shear  forces  for  both  variable 
and  constant  inflow  show  more  than  a  1(X)  per  cent  increase  in  the 
important  third,  fourth,  and  fifth  harmonics  for  the  H-34  at  both  41 
and  112  knots  (Figures  13  and  14).  These  harmonics  are  the  main 
contributors  to  airframe  excitation  due  to  the  filtering  action  of  the 
rotor  hub.  Other  harmonic  shears  are  increased  as  well.  While 
these  increases  in  root  shear  forces  are  substantial,  it  is  estimated 
that  still  larger  forces  are  needed  for  correlation. 


H.  CORRELATION  OF  ONE-HALF  PEAK-TO-PEAK  BENDING 
MOMENTS  ^  _ 


For  the  determination  of  blade  fatigue  life,  knowledge  of  one- 
half  peak-to-peak  vibratory  stresses  is  essential.  The  incorporation 
of  variable  inflow  in  the  analysis  produced  little  change  in  the  ability 
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AIR  LOAD.  LB./m.  AIR  lOAD.  LB./IN. 


FIGURE  9.  H-34  AIR  L0;J)S  AT  /*1  KNOTS 
COMPARISON  OF  CONSTANT  AND  VARIABLE  II^LOW. 
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FIGURE  10.  H-34  FLA1VISE  MOKENTS  AT  41  KIX>TS 
CaiPARISON  OF  CONSTANT  AND  VARIABLE  INFLOW. 
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FIGURE  11.  H-34  AIR  LOADS  AT  112  KNOTS 
COMPARISON  OP  CO^ANT  aND  VARIABLE  INFLON* 
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AZIMUTH  kltaiL,  nvamsa 


FIGURE  12.  H-34  FLATWI2E  MOIfiNTS  AT  112  KNOTS 
COMPARISON  OF  COICSTAI.T  A1;D  VARIABLE  INFLOW. 


HAIU^NIC 

FIGURE  13.  CALCL'UTED  HARMONICS  OF 
ROOT  SHEAR  FORCES  FOR  H-34  AT  41  KNOTS  (CASE  1 ) 
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to  predict  these  stresses.  Typical  results  are  shown  in  the  form  of 
one-haif  peak-to-peak  bending  moment  for  the  H-34  at  41  and  112 
knots  (Figure  15).  At  112  knots  fairly  good  correlation  with  H-34  data 
was  obtained,  but  constant  inflow  gave  somewhat  better  correlation 
than  variable  inflow,  at  least  in  predicting  maximum  values.  At  41 
knots,  the  use  of  variable  inflow  did  improve  correlation.  However, 
the  H-34  data  indicated  higher  values  for  the  moments  than  did  the 
analysis.  It  is  the  high  speed  flight  conditions  which  produce  the  high 
vibratory  stresses  for  which  blades  are  designed  for  fatigue  life.  Con¬ 
stant  inflow  has  produced  good  correlation  of  the  one-half  peak-to-peak 
vibratory  stresses  at  high  speeds,  and  its  continued  use  for  blade  de¬ 
sign  can  be  Justified.  At  high  speeds  the  effects  of  variable  inflow 
are  not  as  significant.  The  main  contribution  to  the  vibratory  stresses 
at  high  speeds  is  made  up  of  the  lower  harmonics.  Constant  inflow 
does  not  provide  the  full  contribution  of  the  higher  harmonics,  but 
does  lead  to  good  definition  of  the  lower  harmonics.  The  addition  of 
the  higher  harmonic  terms  introduced  by  the  use  of  variable  inflow 
may  lead  to  less  accurate  one-half  peak-to-peak  vibratory  stress 
correlation  through  phasing  errors  in  the  summation  of  those  harmon¬ 
ics. 


The  main  contribution  of  variable  inflow  is  in  its  ability  to  pro¬ 
vide  better  definition  of  the  magnitudes  of  the  higher  harmonic  bla^ 
root  shears  which  are  the  forces  producing  airframe  vibration.  These 
forces  are  large  at  low  speeds  as  well  as  high  speed.  The  differences 
in  the  shear  forces  given  by  consunt  and  variable  inflow  are  substan¬ 
tial  as  was  discussed  in  the  previous  section. 

I.  EFFECT  OF  INTERHARMONIC  DAMPING 


The  differences  between  air  loads  calculated  for  the  H-34  at 
70  knots  with  neutral  C.G.  and  at  73  knots  with  aft  C.G.  were  larger 
than  were  expected.  Since  azimuthal  plots  of  the  difference  between 
calculated  air  loads  for  the  two  cases  showed  the  difference  to  be 
primarily  in  the  second  harmonic,  it  was  believed  that  exclusion  of 
interharmonic  damping  might  be  the  cause  of  the  differences,  llie 
reason  is  that  the  main  difference  between  the  neutral  C.G.  case  and 
the  aft  C.  G.  case  is  the  presence  of  first  harmonic  blade  flapping 
about  the  rotor  shaft  axis.  This  flapping  coupled  with  the  cyclic  pitch 
effects  would  cause  an  indicated  difference  in  the  steady  and  secc^ 
harmonic  air  loads.  No  real  difference  should  exist.  It  is  the  role 
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of  the  damping  terms  to  compensate  for  effects  of  flatwise  blade  mo¬ 
tions,  either  due  to  flapping  or  bending,  on  air  loads.  Inclusion  of 
only  the  harmonic  damping  terms  gives  satisfactory  results  for 
neutral  C.G.  cases  where  periodic  flatwise  blade  motions  are  rela¬ 
tively  small.  However,  the  aft  C.G.  case  gives  large  first  harmonic 
flapping  motions  relative  to  the  shaft  axis.  This  creates  substantial 
first  harmonic  damping  forces  and  would  create  large  second  harmon¬ 
ic  damping  forces  if  the  interharmonic  damping  coefficients  were  in¬ 
cluded.  The  important  point  here  is  that  the  aerodynamic  damping 
forces  serve  two  purposes.  First  they  are  corrective  forces  which 
modify  the  undamped  air  loads  calculated  for  a  rigid  non-flapping 
blade  to  give  the  true  air  loads  acting  on  the  flapping  and  bending 
blade  at  its  actual  location  in  space.  These  corrective  forces  do  have 
the  characteristics  of  damping  forces  and,  therefore,  serve  to  damp 
the  motions  which  a  blade  would  have  if  excited  only  by  the  undamped 
rigid- blade  air  loads.  The  two  effects  cannot  be  separated  for  they 
are  in  reality  merely  two  different  ways  of  looking  at  the  same  effect. 
Therefore,  exclusion  of  inierharmonic  damping  forces  can  be  con¬ 
sidered  to  give  either  insufficiently  corrected  air  loads  or  inaccurate 
damping.  This  can  best  be  seen  by  considering  an  example. 

At  a  radial  station  of  r/R  =  .  89,  the  undamped  air  loads  cal¬ 
culated  for  the  neutral  and  aft  C.G.  cases  show  a  difference  (Figure 
16)  which  is  made  up  primarily  of  a  first  and  second  harmonic.  It 
can  be  seen  that  the  first  harmonic  difference  is  virtually  eliminated 
and  only  a  second  harmonic  difference  remains.  The  damping  force 
due  to  first  harmonic  flapping  has  cancelled  the  effect  on  undamped 
air  loads  caused  by  differences  primarily  in  blade  pitch  between  the 
two  cases.  Also  shown  in  Figure  16  is  the  additional  damping  force 
which  inclusion  of  interharmonic  damping  coefficients  for  the  first 
three  harmonics  would  produce.  This  force  would  greatly  diminish 
the  second  harmonic  content  of  the  different  plots  and  make  the  air 
loads  for  the  two  C.G.  cases  match  more  closely,  as  they  should. 

On  this  basis  it  does  apipear  that  the  interharmonic  damping  effects 
are  important  when  aft  C.G.  cases  are  to  be  considered.  Note  that 
if  the  aeroelastic  analysis  had  been  written  in  terms  of  a  tip  path 
plane  axis  system  and  not  a  shaft  axis  system,  the  interharmonic 
damping  effects  due  to  first  harmonic  flapping  would  not  enter  in. 

The  tip  path  plane  would  be  by  definition  the  plane  in  which  no  first 
harmonic  flapping  occurred.  However,  since  the  tip  path  plane  is 
defined  by  blade  response,  it  becomes  an  unwieldy  axis  system  in 
which  to  write  the  equations  of  motion. 
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A  I?.  LOAD 


FIGURE  16.  DIFFERENCES  BETWEEN  CALCUUTED  A  I?.  LO/J)S 
FOR  H-34  at  to  knots,  NEUTRAL  CE?ITni  CF  O'^.'VITr  ;ia) 

w.x't:  ;ft  cf  gravity 
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Imerharmonic  damping  effects  can  be  important  whenever 
there  is  a  large  difference  in  relative  magnitude  between  adjacent 
harmonics  of  blade  motion.  The  first  harmonic  flapping  occuring  for 
the  aft  C.G.  case  considered  here  provides  the  most  obvious  example 
of  such  a  difference,  as  the  second  harmonic  displacement  for  this 
case  is  relatively  small.  Other  instances  of  such  differences  may 
occur  even  for  neutral  C.G.  cases,  for  which  there  is  no  flapping. 
Calculations  of  interharmonic  damping  magnitudes  for  the  first  three 
harmonics  of  the  70-knot  and  73- knot  cases  showed  that  in  the  absence 
of  flapping  the  interharmonic  damping  terms  were  not  erf  sufficient 
magnitude  to  affect  the  airload  signature.  Therefore,  no  significant 
improvement  in  overall  correlation  of  air- load  signature  can  be  ex¬ 
pected  for  neutral  C.G.  cases  as  a  result  of  the  inclusion  of  inter¬ 
harmonic  damping;  however  better  definition  of  some  harmonics  may 
result. 


J.  SUMMARY  OF  DATA  CORRELATION 


The  results  of  this  analytical  study  have  shown  equally  good 
correlation  for  all  air  speeds  considered.  Some  important  differences 
were  noticed  between  neutral  C.G.  and  aft  C.G.  cases.  This  appar¬ 
ently  is  due  largely  to  the  effect  of  interharmonic  damping.  Other¬ 
wise  certain  generalities  can  be  made  about  the  results.  Harmonics 
of  both  air  loads  and  moments  showed  better  amplitude  than  phase 
correlation.  In  fact,  the  azimuthal  plots  show  a  general  trend  to¬ 
wards  a  phase  shift  between  test  and  analytical  data.  This  phase 
shift  is  apparent  in  the  air- load  data  and  is  carried  through  to  the 
moments.  The  similarity  in  correlation  between  test  and  analytical 
air -load  harmonics  and  corresponding  test  and  analytical  bending 
moment  harmonics  points  toward  the  conclusion  that  better  definition 
of  air  loads  is  necessary,  for  bending  moment  correlation  cannot  be 
expected  to  be  better  than  air- load  correlation. 

Another  interesting  result  is  that  a  comparison  between  air 
loads  for  the  H-34  and  the  HU-IA  at  112  knots  indicated  better  air¬ 
load  correlation  for  the  HU-IA.  One  important  difference  between 
the  two  aircraft  is  in  the  number  of  rotor  blades.  The  two-bladed 
HU-IA  could  be  expected  to  show  better  air-load  correlation  as  the 
basis  of  greater  simplicity  of  wake  geometry.  The  four-bladed  H-34 
rotor  may  produce  a  wake  pattern  in  which  interaction  between  rotor 
blades  becomes  more  important.  This  could  lead  to  actual  wake  be- 
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havior  which  deviates  further  from  the  wake  model  used  in  the  analy¬ 
sis  than  does  the  wake  of  a  two-bladed  rotor.  This  observation  is 
merely  a  speculation  at  the  present  time.  However,  it  is  certainly 
obvious  that  the  wake  geometry  considered  in  the  variable  inflow 
program  is  an  idealization  of  true  wake  behavior.  It  is  an  idealiza¬ 
tion  which  makes  possible  a  significant  advance  in  the  ability  to  pre¬ 
dict  air  loads.  Nevertheless,  it  is  not  unreasonable  to  speculate 
that  the  approach  may  be  more  valid  for  a  simpler  two-bladed  rotor 
than  for  a  rotor  with  a  greater  number  of  blades. 


MODIFICATIONS  TO  ANALYSIS 


As  discussed  in  the  previous  section,  further  modifications 
to  the  analysis  are  necessary  to  improve  correlation  with  H-34  test 
data.  Results  indicated  that  better  definition  of  the  air  loads  was 
needed.  Therefore,  three  cases  were  run  to  determine  the  effects 
of  refinements  to  thie  basic  calculation  of  air  loads.  These  are  as 
follows: 


Case  lA  H-34;  steady- state  level  fUght  at  41  knou  using  in¬ 
flow  velocities  generated  from  a  radial  distribution 
of  assumed  downwash. 

Case  4A  H-34;  steady-state  level  flight  at  112  knots  with  air 
loads  recalculated  to  include  effect  of  blade  bending. 

Case  4B  Extension  of  case  4A  to  include  effects  of  blade 
bending  on  inflow  distribution  used  to  define  air 
loads. 


Cases  4A  and  4B  were  run  to  asceruin  the  effect  flexible  mo¬ 
tions  had  on  the  calculation  of  induced  velocities  and  air  loads.  In 
the  aerodynamic  analysis,  the  relation  for  the  vertical  component 
of  air  velocity  on  a  blade  element,  Up,  as  given  belcw,  was  modified 
to  accept  a  radial  and  azimuthal  distribution  of  flatwise  bending  slope, 
p  .  which  was  obtained  from  the  fully  coupled  Myklestad  response 
analysis  for  Case  4. 


Up  =  cos^  (Vi  +  Vsina  t)  Vcosa  isinfi  coe^  (33) 

The  equation  used  in  the  standard  program  is  the  above  equation  with 
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/9  equal  to  the  steady  coning  angle. 

At  the  same  time,  the  angle  of  attack  of  the  blade  at  each 
radial  station  and  azimuth  position  was  modified  by  the  addition  of 
torsional  deformation  of  the  blade  as  obuined  from  the  fully  coupled 
analysis  for  Case  4. 

Air  loads  were  then  recomputed  and  the  trim  conditions  once 
more  satisfied.  This  entire  procedure  can  be  repeated  until  desired 
convergence  is  obtained. 

The  second  modification  at  112  knots  (in  Case  4B)  was  the  in¬ 
clusion  of  blade  deflection  on  variable  inflow  calculations.  The  blade 
displacements  from  Case  4  were  supplied  as  input  to  the  variable  in¬ 
flow  analysis  of  Reference  3.  The  analysis  considers  the  effect  of 
these  motions  only  as  local  effects  along  the  blade.  No  change  in 
wake  behavior  results  from  the  inclusion  of  flexible  motions.  The 
new  values  for  induceu  velocity  which  the  program  calculates  are  then 
used  in  a  second  pass  through  the  Sikorsky  aeroelastic  analysis. 

A  third  change  was  made  for  the  41 -knot  case  (Case  lA)  to 
determine  the  importance  of  wake  variations.  The  wake  transport 
velocity  input  to  the  variable  inflow  program  in  the  other  cases  was 
assumed  to  be  equal  to  the  constant  induced  velocity  determined  from 
momentum  considerations.  It  is  realized  that  this  assumption  is  not 
correct,  for  it  is  known  that  the  rotor  wake  is  definitely  not  comprised 
of  undistorted  helical  vortices  traveling  at  equal  velocities  away  from 
the  rotor.  However,  as  is  stated  in  Reference  3,  the  question  as  to 
what  displacement  time  history  should  be  used  for  a  particular  rotor 
and  operating  condition  remains  unsolved.  Initial  results  presented 
in  Reference  3  also  indicate  that  the  use  of  a  uniform  inflow  distribu¬ 
tion,  determined  from  momentum  considerations,  has  previously  re¬ 
sulted  in  computed  air  loads  which  correlate  well  with  measured  air 
loads  for  a  two-bladed  rotor.  In  addition,  no  provision  is  made  in  the 
existing  variable  inflow  program  for  putting  in  realistic  time  varying 
wake  transport  velocities  which  take  into  account  wake  interaction 
effects. 


The  effect  oi  the  assumed  wake  geometry  on  the  resultant  in¬ 
duced  velocity  and  air- load  distributions  is  of  greatest  significance 
at  low  flight  speed  conditions  because  the  wake  velocity  with  respect 
to  the  rotor,  and  thus  the  distance  between  a  specific  wake  revolution 
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and  the  rotor,  decreases  with  decreasing  flight  speed  and  rotor  angle 
of  attack.  It  was  therefore  reasoned  that  the  results  foi  the  H-34  hel¬ 
icopter  41 -knot  flight  condition  are  the  most  sensitive  to  the  wake 
geometry  assumption. 

An  attempt  had  been  made  to  use  the  variable  induced  velocity 
distribution  from  the  variable  inflow  program  for  the  41 -knot  condi¬ 
tion  as  the  wake  transport  velocity  input  in  a  second  pass  through  the 
program.  This  resulted  in  a  divergent  situation  with  extreme  fluc¬ 
tuations  in  loading  due  to  the  passage  of  blades  close  to  or  actually 
through  vortices  generated  by  previous  blades. 

In  order  to  obtain  additional  information  on  the  effect  of  the 
assumed  wake  geometry  for  the  results  included  herein,  a  radially 
varying  velocity  distribution  was  used  to  define  the  wake  geometry 
(Figure  17).  Tliis  radial  velocity  distribution  consists  of  the  steady 
component  of  the  velocity  distribution  used  in  the  divergent  case  de¬ 
scribed  in  the  previous  paragraph.  The  azimuthal  variation  of  ve¬ 
locity  was  neglected.  With  this  radially  varying  distribution  conver¬ 
gence  was  obtained.  It  is  noted  that  the  average  inflow  value  from 
the  variable  inflow  results  did  not  equal  the  momentum  value.  This 
is  due  to  the  fact  that  the  induced  velocities  from  the  variable  inflow 
program  used  to  define  the  wake  are  the  velocities  that  the  blades 
encounter  and  are  in  the  rotating  coordinate  system,  and  thus  are  not 
the  true  wake  transport  velocities  which  should  be  defined  in  the  n'^- 
rotating  coordinate  system.  To  correct  the  inflow  distribution,  the 
inflow  values  at  each  radial  station  were  decreased  by  three  feet  per 
second  so  that  the  average  value  of  inflow  across  the  rotor  disc 
matched  the  momentum  value. 

The  results  of  the  three  refined  cases  show  chat  ^he  modifi¬ 
cations  produced  little  change  in  either  air  loads  or  bending  moments. 

Case  4A  and  4B  show  only  small  changes  in  results  versus 
azimuth  and  in  magnitude  and  phase  of  harmonic  content. 

Shown  in  Figure  18  is  a  comparison  of  air  load  at  r/R  equal 
to  .  96  for  the  standard  Case  4  versus  Cases  4A  and  4B.  Also  shown 
is  the  result  of  applying  the  variable  inflow  distribution  of  Case  4B 
without  introducing  other  flexibility  effects  on  the  air- load  calculations. 
As  can  be  seen,  small  differences  in  phase  occur  and  the  largest 
change  in  magnitude  occurs  at  around  \f/  =  225  degrees  for  the  com- 
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FIGURt:  17.  RADIAL  DISTRIBUTION  OF  II^CKD  VELOCITY 
USiO)  AS  WAKE  TRAWSFORT  VELOCITY  II^T 
TO  THE  v;.?.:AELr  h-3A 

AT  hi  KI-^5. 
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bined  flexibility  case.  The  fact  that  blade  deflection  has  as  little 
effect  on  air  loads  as  is  shown  is  of  interest.  This  indicates  that 
the  effects  of  higher  harmonic  flapping  on  Up  and  the  effect  of  tor¬ 
sion  on  blade  angle  of  attack  are  only  of  second  order. 

For  Case  lA,  a  significant  change  in  air  loads  occurs  in  the 
180-degree  azimuth,  .  85  blade  station  region.  This  is  due  to  the 
passage  of  the  tip  wake  of  the  previous  blade  directly  beneath  this 
region.  By  assuming  a  greater  wake  velocity  generated  at  the  blade 
tip  the  tip  wake  is  moved  farther  away  from  the  rotor  in  a  given  per¬ 
iod  of  time  and  thus  induces  a  lower  velocity  at  the  blade.  Relatively 
small  changes  are  noted  elsewhere  in  the  rotor  disc,  due  to  the  fact 
that  the  radial  distribution  assumption  introduces  small  differences 
of  induced  velocity  input  in  comparison  to  the  momentum  value  ex¬ 
cept  at  the  tip  as  is  shown  in  Figure  17.  It  is  emphasized  that  this 
simplified  assumption  for  the  wake  transport  velocities  represents 
only  an  initial  effort  to  obtain  some  indication  of  the  influence  of  the 
assumed  wake  geometry  on  the  induced  velocity  distribution  and  as¬ 
sociated  air  loads  and  bending  moments. 

In  general,  the  little  improvement  in  correlation  resulting 
from  these  modifications  indicates  that  more  fundamental  changes  are 
needed,  such  as  the  inclusion  of  wake  contraction  effects  in  the  var¬ 
iable  inflow  analysis. 


EVALUATION  OF  ANALYSIS 


The  results  of  this  program  have  shown  that  fairly  good  cor¬ 
relation  can  be  obtained  from  the  combination  of  variable  inflow  with 
the  fully  coupled  aeroelastic  analysis.  It  is  also  quite  obvious  that 
further  improvement  in  correlation  should  be  possible.  The  process 
of  developing  improvements  for  the  aeroelastic  analysis  and  the  in¬ 
corporation  of  the  variable  inflow  analysis  have  served  to  demonstrate 
the  complexity  of  the  problem.  A  great  number  of  variables  is  in¬ 
volved  in  a  rotor  analysis.  The  success  achieved  with  the  analytical 
methods  chosen  was  in  many  ways  encouraging  and  in  other  ways  dis¬ 
appointing.  In  either  case  much  valuable  information  concerning  the 
importance  of  various  parameters  was  obtained. 

The  basic  question  that  arises  is  what  improvements  should 
be  made  to  improve  correlation.  There  are  three  distinct  areas  to 
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be  considered:  the  basic  blade  dynamic  analysis,  the  trim  conditions 
used,  and  the  air-load  analysis.  The  air-load  analysis  is  the  most 
complex  of  these  areas,  for  it  depends  not  only  on  the  inflow  analy¬ 
sis,  but  also  on  the  way  in  which  the  inflow  is  combined  with  trim 
conditions  and  blade  dynamics. 

The  basic  Myklestad  blade  analysis  used  was  shown  in  SUB¬ 
STANTIATION  OF  METHODS  on  page  3  to  give  valid  results  for  a 
nonrotating  beam.  It  was  found  that  the  number  of  mass  points  cho¬ 
sen  is  adequate  to  define  accurately  the  blade  response.  There  are 
certain  improvements  which  can  be  made  to  this  analysis.  When 
using  an  analysis  based  on  harmonics,  it  must  be  appreciated  that 
the  harmonics  are  not  independent,  but  that  the  use  of  fully  coupled 
motions  introduces  some  interharmonic  coupling  terms.  The  analy¬ 
sis,  as  presently  formulated,  neglects  these  terms  on  the  basis  that 
deflections  are  small,  and  interharmonic  effects  due  to  such  deflec¬ 
tions  can  be  neglected.  It  is  true  that  the  73-lcnot  aft  C.G.  case 
pointed  out  the  importance  of  the  interharmonic  damping  effect  of 
first  harmonic  flapping  when  using  the  shaft-axis  system.  In  that 
case,  the  flapping  motion,  which  is  a  rigid-body  rotation  of  the  blade 
and  not  a  small  deflection,  proved  to  have  an  important  effect.  Con¬ 
sideration  should  be  given  to  extending  the  analysis  to  incorporate 
this  and  any  other  significant  interharmonic  effects.  Effects  to  be 
considered  would  include  such  items  as  the  changes  that  blade  dynam¬ 
ic  rotation  produces  on  flatwise  and  edgewise  bending  through  corre¬ 
sponding  rotation  of  the  principal  axes  of  stiffness  of  the  blade. 

The  basic  trim  procedure  used  is  considered  to  give  a  good 
definition  of  trim  conditions.  Agreement  is  obtained  both  for  gross 
weight,  and  for  the  root -shear  requirements  by  the  double  iteration 
technique  used  for  thrust  moment  and  root  shears.  The  pitch  angles 
obtained  through  use  of  this  procedure  differed  somewhat  from  the 
measured  values  as  shown  in  Table  V  on  page  46.  Collective  pitch 
angle  correlation  showed  the  analytical  values  to  be  generally  low. 
However,  except  for  the  HU-IA  analysis  (see  page  39),  the  radial  dis¬ 
tribution  of  steady  air  loads  correlated  well,  showing  that  the  analy¬ 
tical  thrust  values  matched  test  data.  Cyclic  pitch  values  correlate 
closely,  in  general.  The  differences  which  do  exist  in  cyclic  pitch 
values  and  in  phasing  between  test  and  analytical  blade  response  are 
related.  The  inflow  distribution  as  well  as  rotor  trim  conditions 
affects  pitch  values  and  the  resulting  air  loads.  The  differences  are 
believed  to  be  more  dependent  on  the  inflow  values  than  on  the  trim 


conditions  which  were  used. 


Another  item  which  might  be  considered  as  a  trim  condition 
is  the  determination  of  lag-damper  coefficients.  This  is  most  impor¬ 
tant  for  good  edgewise  moment  correlation.  While  correlation  of 
edgewise  moments  could  be  improved,  correlation  was  found  to  be 
about  the  same  for  the  H-34,  with  lag  dampers,  and  the  HU-IA,  which 
has  no  lag  dampers.  While  this  does  not  verify  the  validity  of  the 
present  lag-damper  coefficient  analysis,  it  does  indicate  that  the  lag 
damper  analysis  is  not  one  of  the  prime  reasons  for  lack  of  correla¬ 
tion. 


The  problem  of  improving  correlation  appears  to  center  on 
the  determination  of  air  loads,  which  is  the  most  complex  part  of  the 
program.  Inaccuracies  in  air-load  calculation  can  occur  from  two 
sources.  First,  the  inflow  velocities  can  be  inaccurate.  Secondly, 
the  way  in  which  rotor  dynamics  is  incorporated  in  the  air  load  cal¬ 
culations  can  introduce  3rrors.  The  differences  which  definition  of 
induced  velocities  can  make  is  quite  obvious  from  comparison  of  the 
results  obtained  using  constant  inflow  with  those  using  variable  in¬ 
flow.  Therefore,  the  degree  to  which  the  variable  inflow  analysis  of 
Reference  3  simulates  the  variable  inflow  distribution  through  a  rotor 
is  most  important.  This  will  be  discussed  in  subsequent  paragraphs. 
Consider  first  the  calculation  of  air  loads  using  a  given  inflow  distri¬ 
bution. 


Initially,  undamped  air  loads  are  calculated  for  the  rotor  sy¬ 
stem  using  the  variable  inflow  distribution  as  well  as  rotor  plane  in¬ 
clinations  and  blade  pitch  determined  from  the  iteration  procedure. 
Airloads  are  corrected  for  blade  flapping  and  bending  motion  effects 
by  the  introduction  of  aerodynamic  damping  terms.  The  importance 
of  including  interharmonic  damping  coefficients  when  using  the  shaft- 
axis  system  has  already  been  discussed.  It  was  shown  that  inclusion 
of  interharmonic  damping  for  the  H-34  73-knot  aft  C.G.  case  would 
make  those  air  loads  compare  with  the  70-knot  neutral  C.G.  air  loads. 
However,  calculated  air  loads  for  70  knots  did  not  correlate  with  test 
data  as  well  as  would  be  desired.  Thus  while  it  would  appear  that 
inclusion  of  interharmonic  damping  would  improve  air- load  definition, 
this  improvement  would  not  be  sufficient  to  achieve  precise  correla¬ 
tion. 


The  effect  of  blade  bending  and  torsion  on  blade  angle  of  attack 
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was  also  introduced  into  the  air -load  calculation  in  the  modified  H-34 
112-knot  case.  However,  the  small  changes  in  results  which  this 
made  would  strongly  indicate  that  this  flexibility  correction  is  not  the 
key  to  better  correlation. 

There  remains  the  complex  area  of  the  variable  inflow  analy¬ 
sis.  Further  improvement  in  the  definition  of  the  variable  inflow  dis¬ 
tribution  over  the  rotor  is  necessary  for  substantial  improvement  in 
correlation.  The  variable  inflow  analysis  of  Reference  3  is  an  ex¬ 
cellent  step  towards  accurate  air- load  definition.  However,  the  model 
used  still  represents  a  substantial  simplification  of  rotor  wake  beha¬ 
vior.  It  is  recognized  that  improvements  in  the  analysis  are  present¬ 
ly  being  investigated.  Such  improvements  should  further  enhance  the 
potential  of  this  analysis. 

It  is  believed  that  full  use  has  been  made  of  the  capabilities 
of  the  present  variable  inflow  analysis  in  this  program.  Further 
manipulations  of  input  values  cannot  be  expected  to  yield  much  im¬ 
provement  in  correlation.  One  of  the  modified  cases,  Case  lA,  used 
a  radial  wake  velocity  variation  with  little  resulting  improvement. 

The  inability  to  obtain  convergence  with  azimuthal  wake  velocity  var¬ 
iation  in  earlier  work  was  unfortunate,  for  it  seems  probable  that 
azimuthal  wake  velocity  variation  would  improve  the  wake  descrip¬ 
tion  and  increase  the  higher  harmonic  content  of  the  air  loads.  This 
area  bears  further  investigation. 

There  are  certainly  numerous  areas  in  the  variable  inflow 
analysis  which  should  be  studied  to  determine  their  importance.  The 
list  of  assumptions  used  in  the  analysis  given  in  THE  CORNELL  VAR¬ 
IABLE  INFLOW  PROGRAM  on  page  27  is  useful  in  this  regard,  for 
each  of  these  assumptions  should  be  evaluated.  Probably  the  area  of 
greatest  concern  is  the  definition  of  wake  geometry.  The  assumption 
that  wake  vortex  elements  retain  the  strength  and  velocity  imparted 
to  them  as  they  leave  the  blade  is  not  realistic.  Wake  vortex  inter¬ 
action  could  be  important.  Wind  tunnel  tests  have  shown  the  impor¬ 
tance  of  hub  and  fuselage  interference  effects  both  on  wake  geometry 
and  on  flow  through  the  rotor.  The  true  wake  geometry  is  more  com¬ 
plex  than  that  which  can  be  derived  with  the  assumptions  used  in  the 
present  analysis.  Further  refinements  do  appear  to  be  necessary  to 
obtain  a  more  realistic  wake  definition. 

The  use  of  vortex  elements  based  on  lifting  line  theory  rather 


than  finite  vortex  and  lifting  surface  theory  may  also  prove  to  be  im¬ 
portant,  especially  in  two  areas.  First,  the  treatment  of  shed  vor¬ 
tices  in  Reference  3  may  not  influence  inflow  correctly  during  the 
important  first  azimuthal  increment  due  to  the  finite  chord  of  the 
blade.  The  use  of  vortex  elements  may  not  give  adequate  definition 
of  the  airflow  behavior  near  the  blade  tip.  The  tip  area  is  most  im¬ 
portant  for  defining  the  response  of  the  entire  blade.  The  sharp  dis¬ 
continuity  at  the  tip  resulting  from  the  lifting  line  theory  may  lead 
to  significant  errors.  The  data  show  that  air- load  correlation  near 
the  tip  does  need  further  improvement. 

To  sum  up,  it  would  appear  that  as  presently  formulated  the 
variable  inflow  analysis  and  the  Sikorsky  aeroelastic  analysis  can 
give  fairly  good  definition  of  rotor  blade  dynamics.  The  improve¬ 
ments  in  correlation  which  were  obtained  make  this  combined  analy¬ 
sis  a  more  useful  design  tool  than  has  been  available.  Both  analyses 
have  room  for  improvement.  However,  significant  improvements 
in  correlation  can  be  made  by  combining  the  present  Sikorsky  aero¬ 
elastic  analysis  with  a  refined  variable  inflow  analysis.  The  develop¬ 
ment  of  the  present  variable  inflow  analysis  by  Cornell  Aeronautical 
Laboratory  has  been  an  important  contribution  to  rotor  system  analy¬ 
sis.  Their  continued  interest  and  the  interest  of  others  in  this  field 
should  lead  to  even  better  definition  of  inflow  in  the  near  future. 
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CONCLUSIONS 


The  purpose  of  this  program  was  to  develop  a  blade  rotor  dy¬ 
namic  analysis  combining  the  Sikorsky  blade  analysis  with  the  CAL 
variable  inflow  analysis  of  Reference  3.  Results  of  this  analysis 
were  correlated  with  flight  test  data.  This  procedure  has  demon¬ 
strated  the  accuracy  which  can  be  obtained  by  this  analysis. 

It  has  also provided  an  opportunity  to  gain  a  better  under¬ 
standing  of  the  significance  of  the  various  parameters  and  methods 
considered.  As  a  result  a  number  of  conclusions  can  be  drawn  from 
this  study.  These  are  summarized  below. 

A.  CORRELATION  WITH  TEST  DATA 

1.  The  combined  Sikorsky-CAL  analysis  provides  an  im¬ 
proved  design  tool  for  the  prediction  of  rotor  dynamic  re¬ 
sponse  for  steady -state  level  flight.  Better  correlation 
with  flight  test  data  was  obtained  using  variable  inflow 
than  has  been  obtained  previously  with  constant  inflow. 

The  improvement  is  most  significant  at  low  airspeeds. 

2.  General  correlation  between  analysis  and  test  was  fairly 
good.  Amplitudes  of  harmonics  of  air  loads,  bending  mo¬ 
ments,  and  torsional  moments  correlated  better  than  do 
azimuthal  phase  angles.  There  is  a  trend  throughout  the 
correlation  study  toward  a  consistent  shift  in  phase  angle 
between  test  data  and  analysis. 

3.  Larger  and  more  realistic  values  for  root -shear  forces 
were  obtained  using  variable  inflow. 

4.  The  variable  inflow  analysis  does  not  match  the  good  corre¬ 
lation  of  peak-to-peak  stress  values  obtained  using  a  con¬ 
stant  inflow  analysis  at  the  high  speed  conditions  used  for 
blade  design. 

5.  Correlation  was  limited  by  the  accuracy  of  air -load  calcu¬ 
lations.  Flatwise  moment  correlation  corresponded  closely 
to  correlation  of  air  loads,  indicating  that  improved  air¬ 
load  definition  would  improve  moment  correlation. 
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6.  Air- load  definition  can  \>e  most  significantly  improved  by 
refinement  in  methods  for  determining  variable  inflow. 


B.  EFFECT  OF  IMPORTANT  PARAMETERS 

1.  The  inclusion  of  blade  bending  deflections  in  the  variable 
inflow  analysis  produces  only  minor  changes  in  inflow  ve¬ 
locities.  These  changes  have  little  effect  on  air  loads. 
Blade  bending  deflections  do  have  an  important  effect  on 
air  loads  through  the  aerodynamic  damping  expressions 

in  the  dynamic  analysis.  The  inclusions  of  blade  torsional 
deformations  and  the  higher  harmonics  of  radial  slope  of 
the  blade  due  to  bending  had  little  effect  on  air  loads. 

2.  Changes  in  first  harmonic  flapping  resulting  from  shifts  in 
C.G.  location  have  little  effect  on  air  loads  for  an  articu¬ 
lated  rotor  with  moderate  offset.  At  present  the  analysis 
does  indicate  a  change  in  air  loads  with  flapping.  This  is 
due  to  the  use  of  a  shaft-axis  reference  system.  The  addi¬ 
tion  of  interharmonic  damping  terms  in  the  analysis  would 
compensate  for  this  change. 

3.  Careful  definition  of  the  physical  and  aerodynamic  repre¬ 
sentation  of  the  blade  tip  area  is  essential.  Closely  spaced 
blade  stations  were  required  near  the  tip  for  adequate 
definition  of  vortex  effects  and  centrifugal  force  effects.  A 
moment  boundary  condition  at  the  blade  tip  was  required  to 
account  for  the  centrifugal  forces  of  the  rubber-  mounted 
blade  counterweights  reacting  against  the  blade  tip  cap  of 
the  H-34.  Such  reactions  are  not  normally  encountered  at 
the  blade  tip. 
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RECOMMENDATIONS 


While  the  results  of  this  program  show  that  improved  correla¬ 
tion  has  been  obtained,  there  obviously  is  still  room  for  improvement. 
The  following  are  some  of  the  areas  which  should  be  considered  for 
further  research. 

1.  Continued  development  of  the  variable  inflow  analysis  to 
provide  refinements  in  inflow  predictions.  Better  defini¬ 
tion  of  wake  geometry  is  an  important  aspect  of  such  work. 

2.  Further  development  of  the  aeroelasdc  analysis  to  consider 
inclusion  of  significant  interharmonic  effects  such  as  inter¬ 
harmonic  damping. 

3.  Incorporation  of  the  combined  analysis  used  in  this  pro¬ 
gram  into  a  fully  integrated  rotor  system  and  airframe 
analysis  using  the  basic  procedure  developed  by  Sikorsky 
Aircraft  (Reference  5).  Rotor  and  airframe  dynamic  re¬ 
sponse,  and  the  interaction  between  the  two  subsystems 
would  then  be  solved  for  simultaneously  by  matching  the 
impedances  of  the  fuselage  hub  and  the  blade  root.  This 
would  provide  a  more  useful  design  tcol  but  would  not  be 
expected  to  provide  improved  air-load  correlation.  The 
shears  at  the  blade  root  however,  should  be  noticeably 
improved. 
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APPENDIX  I 


DEVELOPMENT  OF  CLOSED-FORM  ANALYTICAL  SOLUTION 

- FOR  A  UNIFORM  BEAM - 


The  general  problem  of  a  uniform  pinned- free  nonrotating 
beam  of  mass  ^  per  unit  length  excited  by  a  uniformally  distributed 
p'/rev  harmonic  force,  Fq  Sin  p'u#  t,  is  expressed  by  the  following 
nonhomogeneous  partial  differential  equation. 

a^yrrrr  +  Ytt  "  k'sinp’wt  (34) 

where  a^  =  E  l/fjL  ,  k'  =  Fo//i  .  p'  =  L  2,  ...n  andw  is  the  ex¬ 
citing  frequency.  The  letter  "y"  denotes  the  displacement  of  the  cen¬ 
troid  of  any  section  at  right  angles  to  the  unstrained  central  line. 
Subscripts  r  and  t  indicate  the  respective  differentiation. 

The  equation  rf  flexural  vibration  was  derived  on  the  basis  of 
the  following  assumptions: 

1.  Plane  sections  remain  plane. 

2.  Stress  is  proportional  to  strain. 

3.  Bending  occurs  in  the  principal  plane. 

4.  Slope  of  deflection  curves  are  small. 

5.  Deflection  due  to  shear  is  negligible. 

6.  Rotary  inertia  effect  is  negligible. 

The  above  equation  is  readily  reduced  to  a  homogeneous  or¬ 
dinary  differential  equation  of  amplitude  by  the  steady-state  solution 
in  the  form. 

y(r,  t)  =  X(r)sinp’iit  (35) 

where  X(r)  is  the  spatial  variation  of  the  amplitude. 

Sin  p'ui  t  is  the  steady- state  time  variation  of  the  deflection. 
Substituting  (35)  in  (34),  we  obtain  the  deflection  equation 
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where 


Xrrrr 


X  =  k 


k  =  Fo/ 


anc^  p  =  p'« 


Considering  only  the  homogeneous  portion  of  Equation  (36).  The  com¬ 
plementary  solution  of  the  amplitude  becomes  a  function  of  trigonomet¬ 
ric  and  hyperbolic  functions. 


X(r)c=Asin,^^^  r-f  Bcos^y^/d  r  +  Ccosh ^p/a  r  +  Dsinh ^p/a'  r 

where  A.  B,  C  and  D  are  constants,  whose 
value  depends  upon  the  end  constraints. 


(37) 


The  particular  solution  of  Equation  (36)  can  be  written  in  the  form; 

X(r)  =  ^  k  (38) 

Thus  the  general  solution  of  the  amplitudes  is  defined  as  the 
sum  of  its  complementary  and  particular  form. 


X(r)  =  X(r)c  +  X(r)p 


D  sinh^VpTa  r  -  ^ 


The  boundary  conditions  of  a  pinned-free  nonrotating  beam  are: 

1.  Deflection  y(o,  t)  =  X(o)  sin  p'w  t 

2.  Moment  -El  y,.j.(o.  t)  =  -El  Xj.r(o)  sin  p'«i  t 

3.  Moment  -El  yrr(L.  t)  =  -El  Xj.r(L)  sin  p'«  t 


4.  Shear 


-El  yj-rCL,  t)  =  -El  Xrrr(L)  sin  p'w  t 


rrr^ 


(40) 


Introducing  these  relations  into  (39).  the  constants  are  evaluated 
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A  = 


B= 


^  (a2/p2)  ic 


C= 


D= 


(a2/p2)  k 
1 


{-co8^^^J3  L  sinh^p/a'  L  +  sin^p/a’  L  cosh^p/a'  l1 

(41)  J 
(42) 


(43) 


cosh^ypT?  L  cos L-sinh^p/al-  ^ 


L  sin^p/i  L-sinh^p/a  L  cos.yf/a'"LT 

^  J/ 


(44) 

The  coefficients  being  known,  we  are  able  to  calculate  steady’ 
state  deflections,  and  likewise  by  successive  differentiation  the  steady* 
state  slopes,  momenta  and  shears  produced  by  a  uniformly  distributed 
p'/rev.  harmonic  force. 

By  virtue  of  equation  (35)  the  steady-state  deflection  along  the  nonrota* 
ting  beam  is  given  by: 

y(r.  t)  =  a^  ^  sin^p/S  r  +  ^2  sinh^T^  r  + 

co8.yp/a  r  -  cosh^^p/a  "  2^  sinp’wt 


Differentiating  Equation  (45),  the  slope  becomes: 


(45) 


yv(r.  t)  =  k 


a  co8,^p/a  r  +  **  ’ 

sia^p/a  r  sinh^p/a  r^  sin  p'w  t 
The  spanwise  moment  distribution  is  defined  as: 


(46) 


M  =  E  I  y 


rr 
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M  =  Elak  ^  sin\/p/^'  r  +  sinh\/p/  r  - 

TpH  ^  ^  ^  a 


cos  +  COS  h\/^ «} 

Finally  the  shear  distribution  becomes: 


sin  put 


Q  -  E  I  yrri- 


Q  =  E  1^/7^  k  I"*  j  cosypT^  r  +  <I>2  cosl;^/^  r  + 


r  +  sin 


v^4  sin  p'  u>  t 


(47) 


(48) 


Quantities  and  ^2  appearing  in  Equations  (45),  (46),  (47), 
and  (48)  are  defined  by  the  following  expressions. 


♦.  =  rl  -  s'lVPT^  E  sinlypTJ  L  -  CO  L  cos  V^E-, 


=  rl  -  sir 
\  -  co.^ 


«syb/a  L  siniylj/a'  i.  +  siyt^/a'L  cosh^/a'  E 


(49) 


<X>2  =  f  costiyf)/^  L  co.^yp7^  L  -  L  L  -  1 

I  COS 


'V^i'  ^  *'  ■  ■'  cos,^/f7^ 


(50) 
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FIGURE  25.  CASE  1  :  H-34  CHORDWISE  *40MENTS  AT  41  KNOTS 
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FIGURE  35.  CASE  i*:  H-3A  CHORDWISE  MOMENTS  AT  112  KNOTS 
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